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1.0  INTRODUCTION 


1.1  BACKGROUND 

The  excitation  of  gaseous  atomic  and  molecular  species  by  high- 
energy  electron  impact  has  long  been  of  fundamental  interest.  However, 
during  recent  years  additional  attention  has  been  given  this  subject  be¬ 
cause  of  the  applications  of  such  knowledge  including  the  use  of  electron 
impact  excitation  for  determining  various  properties  of  the  gaseous 
species.  Our  own  interest  is  motivated  by  the  use  of  electron  beam 
fluorescence  diagnostics  of  high-speed  flow  fields  of  gases  and  their 
mixtures  for  the  purpose  of  determining  local  values  of  gas  specie 
density  and  temperature  by  observing  the  fluorescence  radiation  result¬ 
ing  from  electron-molecule  collisions  of  5  to  50  keV  energy.  Although 
the  breadth  of  the  energy  range  of  interest  is  somewhat  arbitrary, 
various  application  restraints  do  exist  that  place  practical  limits  on 
both  the  lower  and  upper  bound  of  the  electron  energy.  The  primary 
constraint  encountered  in  diagnostic  applications  is  that  of  an  electron 
beam  spreading  along  the  injection  path  as  a  result  of  both  elastic  and 
inelastic  scattering  processes,  thereby  degrading  the  spatial  resolution. 
Consequently,  on  this  basis,  one  would  desire  to  employ  beam  energies 
on  the  order  of  50  keV  to  minimize  spreading  problems;  energies  greater 
than  about  50  keV  are  obviously  practical,  but  little  is  gained  in  the  area 
of  beam  spreading  by  so  doing,  and  such  high  energies  are  not  deemed 
worth  the  added  arcing  and  corona  discharge  problems  encountered. 

A  second,  but  less  important,  consideration  is  encountered  when 
applying  this  technique  to  the  measurement  of  low  gas  density  samples 
as  are  typically  found  in  the  far  field  of  hypersonic  expansions.  For 
these  cases,  one  desires  as  much  radiative  intensity  as  is  possible 
within  the  spatial  resolution  requirements,  and  consequently,  the  use  of 
lower  beam  energies  and,  thereby,  larger  cross  sections  is  attractive. 

A  lower  limit  of  1  to  5  keV  for  electron  beam  energy  in  these  instances 
is  not  unreasonable. 

As  will  be  discussed  in  detail  in  a  following  section,  the  electron- 
molecule  excitation  process  can  be  categorized  according  to  the  changes 
in  angular  momentum  and  spin  quantum  numbers  resulting  from  the 
collision.  Generally,  such  collisions  are  denoted  as  allowed  or  for¬ 
bidden,  and  the  degree  of  forbiddenness  can  be  further  refined  as  for¬ 
bidden  in  orbital  angular  momentum  or  spin  angular  momentum  or  both, 
and  each  type  of  transition  is  experimentally  characterized  by  its 
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dependence  on  electron  energy.  As  is  well  known,  in  general,  the 
allowed  transitions  are  much  stronger  than  those  which  are  first-order 
forbidden,  and  as  a  result,  allowed  transitions  are  usually  employed 
for  diagnostics  applications.  However,  even  though  allowed  transitions 
are  preferable  for  use,  one  cannot  ignore  the  weaker,  forbidden  tran¬ 
sitions  for  many  applications.  As  an  example,  if  one  is  studying  the 
flow  field  of  a  binary  mixture  of  He  and  N2  for  which  the  He  mole  frac¬ 
tion  is  small,  the  strong,  allowed  He  transition  S1?-*  21S  at  5016  A 
can  easily  be  obscured  by  the  Nil  transition  3F°->  3D  at  5016.  4  A, 
which  is  the  result  of  the  dissociative -excitation  of  N2  by  electron  im¬ 
pact.  Obviously,  both  for  experiment  design  as  well  as  data  interpre¬ 
tation  the  dependence  of  the  Nil  line  excitation  cross  section  with  elec¬ 
tron  energy  as  well  as  its  absolute  value,  density  and  current  depend¬ 
ences  are  desired.  A  related  question  exists  when  one  considers  flow 
fields  containing  the  molecule  and  at  least  one  of  its  constituent  atomic 
species,  such  as  N2/N,  NO/N  or  O,  H2/H,  HF/H  or  F,  and  F2/F  mix¬ 
tures.  Obviously,  to  determine  the  density  of  the  atomic  specie,  say 
N,  radiative  transitions  of  the  NI  or  Nil  species  must  be  observed; 
higher  ionized  species  such  as  NIII,  NIV,  etc.  ,  are  quite  unlikely. 

Since  the  electron  impact  on  both  N  as  well  as  N2  produces  NI  and  Nil 
radiation,  one  desires  to  know  to  what  extent  the  molecular  and  atomic 
species  each  contribute  to  the  resulting  composite  radiation.  To  do  so, 
one  must  know  the  cross  section  for  the  production  of  each  specie  as 
well  as  its  energy  dependence.  For  the  atomic  specie,  inspection  and 
comparison  of  the  ground  state  quantum  number  description  with  that 
of  the  excited,  radiating  state  can,  in  many  cases,  give  one  the  gross 
characterization  of  the  excitation  process  as  either  allowed  or  forbidden. 
The  situation  is  less  simple  for  dissociation  processes  for  molecules, 
since  the  process  involved  is  at  least  three-body  in  nature. 

This  report,  which  is  the  first  in  a  series  of  studies  of  dissociative- 
excitation  processes  of  electron-molecule  collisions,  describes  the  re¬ 
sults  of  the  measurement  of  the  energy  dependence  of  the  dissociative - 
ionization  cross  sections  for  the  excitation  of  the  Nil  states  3F°  and 
by  electron  impact  on  both  N2  and  NO.  The  results  were  anticipated  to 
characterize  each  transition  as  to  whether  it  is  allowed  or  forbidden  and 
to  determine  the  degree  of  linearity  of  the  radiative  intensity  with  gas 
density  and  electron  beam  current.  Finally,  the  study  was  expected  to 
show  to  what  extent,  if  any,  the  excitation  cross  sections  were  affected 
by  chemical  binding.  As  the  following  sections  show,  these  goals  were 
achieved. 
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1.2  THEORY 

The  three  types  of  electron  excitation  processes  of  interest  for  this 
work  are  as  follows: 

Molecular  ionization. 


+  N 2  -*  (N|)»  +  2e~ 

L  N$  +  hu  (1) 

Homonuclear  molecular  dissociative  ionization, 

e"  +  N2  -  (N+)**  +  N  +  2e- 

L  (N+)*  +  hi/  w 

Heteronuclear  molecular  dissociative  ionization 

e-  +  NO  -*  (N+)**  +  0  +  2e~ 

L  (N+)“  +  hi/  (3) 

where  (  )*"■  denotes  an  excited  state  of  the  parenthetically  designated 
specie. 

For  the  processes  represented  by  Eqs.  (1)  to  (3),  it  is  assumed 
and  the  experiment  is  configured  such  that  no  collisional  quenching 
effects  exist. 

The  radiative  transitions  of  interest  are 

n2  b2s;  -  n2  x2s; 

for  the  process  of  Eq.  (1)  and 


NIl’D  ■»  NI^P0 


and 


NII3F°  -  NII3D 
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for  the  processes  of  Eqs.  (2)  and  (3).  The  following  paragraphs  in  this 
section  describe  the  theoretical  basis  for  the  description  of  electron/ atom- 
molecule  collisions,  and  the  section  concludes  with  the  application  of  the 
theory  to  the  processes  of  Eqs.  (1),  (2),  and  (3). 

1.2.1  Differential  Cross  Section 


The  present  derivation  and  discussion  essentially  follows  that  of 
Inokuti  (Ref.  1).  Consider  an  electron  of  velocity  v  which  collides  with 
a  stationary  atom  in  its  ground  state  and  is  inelastically  scattered  in 
the  spherical  polar  direction  (9,  <f>)  of  the  center  of  mass  system  into  the 
solid  angle  element  (du),  and  assume  that  the  atom  undergoes  a  tran¬ 
sition  from  the  ground  state  (o)  to  state  (n)  of  excitation  energy  (En). 

When  the  incident  electron  is  sufficiently  fast,  the  collision  may  be 
regarded  as  a  sudden  and  small  external  perturbation  on  the  target  atom. 
The  differential  cross  section  (do^),  calculated  in  the  first  Born  approxi¬ 
mation,  is  then  given  by  (Ref.  1): 


dan  =  (M2/4ff?fi4)(k  vk)  |f  exp  (iK  •  ?)<£*(?,  .  .  .?z)V0o(i 


r7)dZr  dr 


Yd co 


(4) 


where  M  is  the  reduced  mass,  r  is  the  position  of  the  incident  electron 
relative  to  the  center  of  the  atom,  -fik  and -hie  are  the  initial  and  final 
momenta  of  the  incident  electron,  'hK  =  h(k  -  k'  )  is  the  momentum 
transfer,  and  the  t//‘s  are  the  atomic  eigenfunctions  as  functions  of  the 
coordinates  rj  of  the  atomic  electrons.  The  3Z-dimensional  volume 
element  is  denoted  as 


dZr  =  drj  dr2  di\j  ■  •  .  d?z 

The  interaction  energy  V  is  the  Coulomb  interaction  between  the 
incident  and  atomic  electrons. 


V 


(5) 


The  nuclear  interaction  energy  has  been  omitted  in  Eq.  (5)  since  it 
yields  no  contribution  to  Eq.  (4)  due  to  the  orthogonality  of  initial  and 
final  states. 
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By  using  the  Bethe  relation. 


-l 


/  exp  (iK  •  r)  | 


r  - 


d?  =  477  K~2  exp  (i  K 


(6) 


Eq.  (4)  becomes 


d an  =  (4M2e4^4K4)(k'/k)ifa(K)|2d<u 


(7) 


where  en(K)  is  the  atomic  matrix  element 


fn(K)  =  <n 


2  exp  (iK  ■  r.) 
j=l  J 


0>  =  2  exp  (iK  •  r.)^D  d' 

i=l  ‘ 


(8) 


Generally,  one  may  consider  dan  as  independent  of  <j>  because  the 
ground  state  is  spatially  symmetric  or  because  the  atoms  are  oriented 
at  random  and  an  average  of  all  atomic  orientations  is  implied.  Thus, 
|en(K)|2  is  a  function  of  a  scalar  variable  K,  and  one  writes  |en(K)|2 
instead  of  |en(K)j2.  Further,  when  one  expresses  dcrn  in  terms  of  dK, 
one  implies  integration  over  <j>  since  K  is  independent  of  <j>.  Therefore, 
one  replaces  du  by  2ir  sin  9d6  =  w(kk')"^  d(K2)  in  Eq.  (7)  to  obtain 


dcr  =  ('W M2c4/'fl4K4k2)  )e  (K),2  d(K2) 


(9) 


The  quantity  en(K)  reflects  the  dynamics  of  the  atom  and  is  known 
as  the  inelastic-scattering  form  factor.  For  our  purposes,  however, 
one  more  often  uses  the  generalized  oscillator  strength: 


fn(K)  =  (2mEn/*2K2)  |,b(K)| 


(10) 


where  m  is  the  electron  mass.  By  using  the  Bohr  radius 
aQ  =-fi2/mc2  =  0.  52918  x  10"®  cm  and  the  Rydberg  energy 
R  =  me4/2*2  =  13.  606  eV,  Eq.  (10)  becomes 


'.<»  -  (U) 

The  generalized  oscillator  strength  is  a  straightforward  generalization 
of  the  optical  (dipole)  oscillator  strength  fn,  defined  as 

-  «Vb>m„2  (12) 
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where 


M2 

n 


n 


/ 


V 

j=l 


Xj^0dZr 


(13) 


is  the  dipole -matrix -element  squared  and  xj  is  a  component  of  rj.  The 
relationship  between  the  dipole -matrix  elements  and  types  of  transitions 
is  discussed  at  the  end  of  this  section.  By  expanding  the  exponential  of 
Eq.  (8)  into  the  familiar  power  series  and  using  the  orthogonality  of  the 
initial  and  final  atomic  eigenstates,  it  can  be  shown  that 


lim  f  (K.)  =  f 

„  n  n 

k.-*o 


(14) 


Equation  (14)  is  especially  important  because  it  relates  the  collision  of 
fast  charged  particles  with  photoabsorption  processes. 

When  the  incident  particle  is  an  electron,  mv^/2  =  E  represents 
the  kinetic  energy,  and  Eq.  (9)  becomes 


4ff a”  f  (K) 

-  Ti  f-R  d[£n(Kao)  1 


(15) 


Noting  that  M/m  =  1  applies  approximately  for  incident  electrons, 
and  using  the  energy  conservation  relation. 


<*k)2  =  (fit') 2  +  2MEn 

and  the  definition  of  K  =  k  -  k",  we  find  that 


(16) 


(17) 


For  a  given  value  of  En,  (Ka0)^  has  a  minimum  at  6  =  0  and  a 
maximum  at  0  =  it.  For  En/E  «  1  (Born  approximation),  one  obtains 
from  Eq.  (17) 


<k>L„ 


<*•->««  =  4 


4  [>  •  j©  -  •(©■)] 


(18) 

(19) 
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For  excitation  to  continua  (e.  g. ,  ionization),  the  excitation  energy 
is  no  longer  a  discrete  variable  but  is  a  continuous  variable  E'  taking 
all  real  values  greater  than  the  first  ionization  threshold  Ij.  Now,  the 
cross  section  (at  a  fixed  E)  for  excitation  to  continuum  states  between 
E'  and  E'  +  dE'  is  (dar/dE')dE',  and  dcr/dE'  is  the  density  of  the  cross 
section  per  unit  range  of  E'.  Alternatively,  a(E')  may  be  considered 
as  the  cross  section  for  excitation  to  all  states  up  to  E'  and  dcr/dE'  as 
the  derivative  of  cr(E').  Thus  Eq.  (15)  may  be  adapted  to  give  the  differ¬ 
ential  of  dcr/dE'  as 


E/R 


R 

E' 


df(K.E') 

dE' 


d[MKa0)2] 


(20) 


The  density  of  the  generalized  oscillator  strength  per  unit  range  of  E', 
df(K,  E')/dE',  is  defined  by 

df(K,E')/dE'  =  2(En'R)r|fn(K)|2/(Kao)2ld(En  -  E') 

n  L  J  (21 


where  the  summation  runs  over  all  excited  states— discrete  as  well  as 
continuum. 


It  is  well  known  from  quantum  mechanical  perturbation  theory  that 
the  probability  for  a  transition  from  a  state  m  to  a  state  n  is  proportional 
to  the  factor 


<n  |e 


ik 


P  1  m>  | 2  =  1  f  e 


ik  ■ 


rfi 


P^n 


dr 


(22) 


where  the  tp's  are  the  atomic  eigenfunctions,  k  is  the  wave  propagation 
vector,  and  n  •  p  is  the  component  of  the  electron  momentum  p  along 
the  direction  of  polarization  of  the  incident  radiation,  i.  e. ,  the  direction 
of  the  electric  field  vector.  Usually  the  wavelength  of  the  emitted  or 
absorbed  radiation  is  much  greater  than  the  dimensions  of  the  atom, 
whereas  the  wave  functions  are  nonvanishing  only  in  a  region  of  the  size 
of  the  atom.  Thus,  the  quantity  k  •  r  in  Eq.  (22)  is  small  compared  to 
unity  and  from  the  expansion 


.ik 


•  r  =  1  -  i(k  .  7) 


[-Kk  .7  )Y 
2< 


(23) 
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-*  ^ 

we  can  replace  exp  <ik  •  r )  by  1.  The  probability  for  a  transition  from 
a  state  m  to  a  state  n  is  then  approximately  proportional  to 

Pnm  *  M  "  ‘  Plm>|2  =  |n  •  <n|p|m>|2  (24) 

If  the  matrix  element  in  Eq.  (24)  vanishes,  the  second  term  in  the  ex¬ 
pansion  (Eq.  (23))  must  be  considered.  By  using  the  unperturbed 
Hamiltonian  and  the  commutation  relations  for  momentum  and  coordi¬ 
nate  operators,  it  can  be  shown  in  general  that  the  matrix  element  of 
the  momentum  operator  in  Eq.  (24)  can  be  replaced  by  the  matrix  ele¬ 
ment  of  the  coordinate.  As  the  dipole  moment  is  usually  defined  by 
e  <n|r|m>,  it  follows  that  the  probability  for  a  transition  from  a  state 
m  to  a  state  n  is  proportional  to  the  dipole -matrix- element  squared, 
i.  e. , 


Pnm  K  I"  1  <nlr|m>l  (25) 

For  a  multi-electron  atom,  Eq.  (25)  becomes 

Pnm  *  l<n|M|m>|2  (26) 

where  M  is  a  vector  with  components 

jJ'i’V  Jl'l’"!- 

and  the  atomic  eigenfunctions  are  denoted  by  |n>  =  IS>IL>ln>  and 
|m>  =  ISG>  tL0>ln0>  for  L-S  coupling  and  by  |n>  =  |J>  |n>  and  tm>  =  Uo>|n0> 
for  J- J  coupling.  The  total  quantum  numbers  (n,  nG)  identify  the  final 
and  initial  electronic  states,  respectively.  The  existance  of  nonvanish- 
ing  values  of  the  dipole -matrix- element  of  Eq.  (26)  determines  the 
selection  rules  for  electric  dipole  radiation.  If  either  type  of  coupling 
is  assumed,  the  selection  rule  for  the  total  angular  momentum  is 
AJ  =  0,  ±1  with  the  restriction  that  J  =  0-*J=0is  not  permitted.  In 
addition,  the  initial  and  final  states  must  have  opposite  parity.  For 
L-S  coupling,  there  is  a  strong  interaction  between  individual  electronic 
spin  and  orbital  angular  momentum  vectors  to  give  the  resultant  L  and  S 
with  the  selection  rules  AS  =  0,  and  AL  =  0,  ±1.  For  J-J  coupling,  which 
is  generally  prevalent  for  large  Z  species,  there  is  a  strong  interaction 
between  the  spin  and  orbital  angular  momenta  of  the  individual  electron 
giving  the  resultant  total  electronic  angular  momentum  J  with  the  selec¬ 
tion  rule  AJ  =  0,  ±1  for  the  electron  whose  state  is  changed. 
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Transitions  which  obey  the  selection  rules  described  above  are 
called  electric  dipole  or  optically  allowed  transitions,  whereas  those 
transitions  which  violate  these  rules  are  optically  forbidden.  If  both 
the  initial  and  final  state  wavefunction  are  spherically  symmetric,  the 
integral  in  Eq.  (22)  is  identically  zero,  and  transitions  between  these 
states  are  strictly  forbidden.  The  selection  rules  for  forbidden  tran¬ 
sitions  can  be  determined  by  evaluating  the  matrix  elements  obtained 
by  including  the  second  term  of  the  expansion  Eq.  (23).  Note  that  the 
selection  rules  presented  here  are  all  first  order  in  that  the  probability 
for  the  transition  was  obtained  from  first-order  perturbation  theory. 

1.2.2  Integrated  Cross  Section 

The  cross  section  for  excitation  to  a  specific  state  n  of  an  atom  or 
molecule,  regardless  of  the  angle  of  scattering  of  an  incident  particle, 
is  defined  as  the  integrated  cross  section  of  the  transition.  For  the 
nonrelativistic  case,  crn  is  simply  an  integral  of  dcrn,  as  given  by 
Eq.  (15),  over  all  kinematically  possible  values  of  the  momentum  trans- 
fer -ttK,  i.  e. , 


a 


n 


4/73  2 
o 

E/R 


/ 


o  nun 


f„(K)  d<Kao)2 
T, 7 «  (Kb/ 


(27) 


2  2 

where  (Kaolin  and  (Ka0)max  are  given  by  Eq.  (18)  and  Eq.  (19),  re¬ 
spectively. 


Since  the  Born  approximation  assumes  sufficiently  large  incident 
electron  kinetic  energy,  Bethe  recognized  the  usefulness  of  expressing 
the  cross  section  in  terms  of  an  asymptotic  expansion  in  inverse  powers 
of  electron  kinetic  energy.  To  explain  the  general  features  of  the 
asymptotic  expansion,  we  consider  the  generalized  oscillator  strengths 
of  the  optically  allowed  helium  2*P  excitation  (Ref.  2).  A  plot  of 
fn(K)/(En/R)  versus  in  (Ka0)2  is  given  in  Fig.  1.  The  integral  in 
Eq.  (27)  is  represented  by  the  area  under  the  curve  limited  by 
in  (Kao)2^  and  in  (KaoJ^iax*  Because  the  integrand  decreases 
rapidly  for  large  values  of  (Ka0)2,  for  large  values  of  E,  the  upper 
limit  of  the  integration  may  be  accurately  taken  as  being  infinite. 

Further,  (Kao)^^  in  in  general  small  and  decreases  with  increasing  E. 
Now  choose  the  value  (Ka0)2  such  that  the  areas  ABC  and  CDF  in  Fig.  1 
are  equal.  Thus  it  becomes  obvious  that  the  required  area  is  given  in 
the  Born  approximation  (large  E)  by  the  product  of  in  (Ka0)2  -  in  (Kao)2^ 


15 


AEDC-TR-75-77 


and  the  dipole -matrix- element  squared  =  fn  R/En  which  is  defined 
by  Eqs.  (12)  and  (13).  Analytically, 

tn(fe0>2  =  C  d[Fn(Kaol2]  -  f°  |l  -  ^1  dtFn(Ka0)2] 

°  o  fn  L  n  J  (28 

—  2 

Thus,  (Ka0)  is  independent  of  E,  and  Eq.  (27)  becomes 

EnlKaJ® 

n  >•  o  m  a  x 

W  I  fn<K)  O 

•.-■otJ  ,  yji  ^ 

Pn(Ka  )‘in 

o  m  id 


ofe)j 

)  "  (Ka)2. 

I  o'  min  1 


(Ka0)2 


Figure  1.  Momentum  dependence  of  the  generalized  oscillator  strength 
for  the  optically  allowed  helium  21  P  excitation. 
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Following  the  practice  of  other  authors  in  recent  years,  we  define 


£nCn  = 


(30) 


Consequently,  the  Bethe  asymptotic  formula  for  optically  allowed  ex¬ 
citation  by  electron  impact  is  given  by 


47ra2n 

O 


K  «CnE> 


(31) 


The  optically  forbidden  excitation,  for  which  fn  =  0,  follows  in  an 
obvious  manner.  A  plot  of  fn(K)/(En/R)  versus  jgn(Ka0)2  for  the  helium 
(2*S)  excitation  is  shown  in  Fig.  2.  The  asymptotic  cross  section  for 
optically  forbidden  excitations  is  simply 

a2R 

an  =  ~T~  b«  (32) 


where  bn  is  defined  by 


—  o© 


fn(K) 

F.n/R 


d[fn(Kao)2] 


(33) 


Note  that  the  constants  bn  and  Cn  are  each  properties  of  the  respective 
transition  to  state  n  and  depend  upon  the  excitation  energy  En. 


Figure  2.  Momentum  dependence  of  the  generalized  oscillator  strength 
for  the  optically  forbidden  helium  21  S  excitation. 
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For  excitation  to  continua,  the  integrated  cross  section  do/dE'  per 
unit  range  of  excitation  energy  E',  using  the  differential  generalized 
oscillator  strength  df(K,E')/dE'  in  place  of  fn(K),  is  defined  by  an  equa¬ 
tion  similar  to  Eq.  (27).  The  Bethe  procedure  then  leads  to 

fe^)(#)(dS)«cE-E)  (34) 

where  df/dE'  =  [df(K, E')/dE'  Jj^=q  is  the  differential  optical  oscillator 
strength  and  Cg'  is  defined  by  Eqs.  (28)  and  (30)  after  replacement  of 
fn(K)  by  df(K,E')/dE'. 

The  ionization  cross  section  for  an  atom,  which  is  obtained  by  inte¬ 
grating  Eq.  (34)  over  the  continuum  energy  E',  is  given  by 


where 


(in  r\ 

«l  -  i~)  M?  e"(CiE) 


(35) 


(36) 


and 


*1  (37) 

The  quantity  M^,  which  is  related  to  the  distribution  of  the  optical 
oscillator  strength,  is  called  the  effective  dipole -matrix-element 
squared  for  ionization.  Since  various  processes' can  occur  as  a  result 
of  energy  transfer  E'  in  a  single  collision,  the  quantities  M2  and  C*  are 
expressed  containing  n(E'),  which  is  defined  as  the  efficiency  for  ioni¬ 
zation  at  excitation  energy  E'. 


1.2.3  Molecular  Excitation  Cross  Sections 


When  applying  the  differential  cross  section  formula  Eq.  (15)  to 
molecules,  the  symbol  o  or  n  represents  a  set  of  quantum  numbers 
(electronic,  vibrational,  and  rotational)  which  designates  a  molecular 
state,  and  the  quantity  |en(K)l2  requires  further  specification  to  include 
the  molecular  internal  degrees  of  freedom.  For  a  diatomic  molecule. 
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the  final  state  is  specified  by  the  electronic  quantum  number  n,  vib¬ 
rational  quantum  number  v,  and  rotational  quantum  numbers  J  and  M, 
with  corresponding  initial  state  quantum  numbers  n0,  vQ,  Jo*  and  Mo. 
By  assuming  the  Born-Oppenheimer  separation  between  the  electronic 
and  nuclear  motions  and  neglecting  the  rotational-vibrational  coupling, 
the  final  state  wave  function  is 


and  that  of  the  initial  state  is 


O  O  *0  0 


(38) 


(39) 


where  the  fs  are  electronic  wavefunctions  depending  on  the  coordinate 
?i(i  =  1,  2,  ...  Z)  of  all  the  molecular  electrons  and  on  the  nuclear 
coordinates  p(p,  0,  o),  the  X’s  are  the  vibrational  wavefunctions  depend¬ 
ing  on  the  internuclear  distance  p,  and  the  Y’s  are  spherical  harmonics. 
The  form  factor  for  the  transition  is  then  given  by 


f(K;n,v,J,M  -  n0,v0,J0.M0)  =  /X^<P>Y;M(8.«£n(K;p,0,*)  ^^(©.‘VdptKcos  9)dO 


(40) 


where  en(K;p,  q,  <j>)  is  the  electronic  part  defined  by 

z 


€n(K.;p,0,<J>)  =  S  exp  (iK  •  7.)«fn 

j=  1  J  0 


(41) 


Since  rotational  level  spacings  of  molecules  are  much  smaller  than  the 
energy  resolution  of  most  collision  experiments,  they  are  treated  as 
effectively  degenerate  and,  therefore,  summed  over  all  J  and  M  and 
averaged  over  JQ  and  M0.  The  resulting  squared  form  factor  for  the 
transition  from  the  vibrational  level  v0  in  the  initial  electronic  state  nQ 
to  the  vibration  level  v  in  the  final  electronic  state  n  is  then  obtained  by 
using  the  closure  property  of  spherical  harmonies  as 

|«n(K;v  -  v0)|2  =  (4n)"1  JJd(cos  0)d$  |/X*  v(p)en(K;p,0,O)  Xn^v{p)p2ip\ 2  (42) 

o 

Note  this  result  shows  that  |en(K;v  v0)l  is  equal  to  the  average  of  the 
squared  form  factor  evaluated  for  all  possible  molecular  orientations. 
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Since  en(K;p,  g,$)  is  generally  a  more  slowly  varying  function  of  p 
than  are  the  vibrational  wave  functions,  en(K;p,  9,  $)  may  be  replaced 
by  its  value  en(K;/3,  ©,  <t>)  at  a  representative  p  value  of  p,  which  may  be 
taken  as  the  equilibrium  internu clear  distance  pe  if  vQ  =  0.  Physically, 
this  means  that  the  nuclei  do  not  move  appreciably  in  the  period  during 
which  electrons  undergo  a  transition  under  the  influence  of  the  incident 
electron.  Thus,  Eq.  (42)  becomes 

|«b(K;v  -  v0)|2  =  (-Itt)'1  /fd(cos  0)d4>  |fnK;p,0,<t>)|  2  |  fX*v(p)X„^v^(p)p2dp|  2  (43) 

Note  that  the  first  factor  in  Eq.  (43)  deals  primarily  with  the  electronic 
motion,  and  thus  may  be  called  the  electronic  form  factor  squared; 
whereas  the  second  deals  exclusively  with  the  nuclear  vibration  and  is 
identical  to  the  Franck -Condon  factor  and  is  well  known  in  molecular 
spectroscopy.  Applications  of  the  Franck- Condon  principle  will  be 
shown  later  in  this  section. 

The  differential  cross  section  equation  for  molecular  electron  ex¬ 
citation  to  state  n  is  the  same  as  in  the  atomic  case  (Eq.  (15))  with  the 
form  factor  in  the  generalized  oscillator  strength  (Eq.  (11))  redefined 
by  Eq.  (43)  and  n  redefined  to  include  the  entire  set  of  molecular  quan¬ 
tum  numbers.  Applying  this  reinterpretation  of  the  form  factor  and  n 
to  Eq.  (20)  gives  the  differential  of  the  cross-section  density  per  unit 
range  of  E'  for  excitation  to  continua. 

The  integrated  cross  section  for  molecular  excitation  to  state  n, 
where  n  contains  the  required  set  of  molecular  quantum  numbers,  is 
given  by  Eq.  (27)  if  the  generalized  oscillator  strength,  or  the  form 
factor,  is  redefined  by  Eq.  (43).  The  Bethe  procedure  then  leads  to 
an  equation  identical  in  form  to  Eq.  (31)  for  optically  allowed  transitions 
and  to  Eq.  (32)  for  optically  forbidden  transitions.  The  selection  rules 
for  these  transitions  are  obtained  just  as  in  the  atomic  case  by  evaluating 
the  matrix  elements  of  the  electric  dipole  moment  and  are  described  in 
detail  by  Herzberg  (Ref.  3). 

For  excitation  to  continua,  the  integrated  cross  section  do/dE'  per 
unit  range  of  excitation  energy  E',  using  the  reinterpreted  differential 
generalized  oscillator  strength  for  molecular  processes  df(K,E')/dE' 
instead  of  fn(K),  is  defined  by  an  equation  similar  to  Eq.  (27).  The 
Bethe  procedure  then  leads  to  an  equation  identical  in  form  to  Eq.  (34) 
for  excitation  to  continua. 

For  an  atom,  excitation  is  generally  to  a  specific  state  by  an 
allowed  or  forbidden  transition  or  to  continua  as  an  ionization  process 
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with  integrated  cross  sections  described  by  Eq.  (31),  (32),  and  (35), 
respectively.  For  a  molecule,  additional  excitation  channels  are  open 
for  each  electron-molecule  collision;  these  can  be  summarized  by  three 
basic  diagrams  that  also  serve  to  illustrate  the  consequence  of  the 
Franck-Condon  principle  (Ref.  3).  In  Fig.  3,  the  lower  curves  repre¬ 
sent  the  initial  state  of  the  molecule.  By  assuming  the  Franck-Condon 
principle,  the  transition  takes  place  so  quickly  that  the  nuclei  do  not 
have  time  to  move  an  appreciable  distance,  and  since  the  nuclear  sepa¬ 
ration  in  the  lower  vibrational  level  will  effectively  lie  between  a  and  b, 
the  upper  state  of  the  molecule  will  be  represented  by  points  lying  be¬ 
tween  c  and  d  on  the  upper  curve. 


Figure  3.  Molecular  electronic  transitions  illustrating  the 
consequences  of  the  Franck-Condon  principle. 

Figure  3a  illustrates  excitation  to  a  stable  molecular  state,  and 
Fig.  3c  shows  excitation  to  continua  resulting  in  dissociation  of  the 
molecule,  while  Fig.  3b  is  a  mixture  of  excitations  leading  to  stable 
molecules  and  to  dissociation.  Similarly,  ionization  of  a  molecule 
can  be  described  using  Fig.  3  where  the  upper  potential  energy  curve 
corresponds  to  an  electronic  state  of  the  molecular  ion.  When  atomic 
excitations  resulting  from  the  dissociation  of  a  molecule  are  of  interest, 
many  more  excitation  channels  become  available  thus  further  compli¬ 
cating  the  analysis. 

Fortunately,  these  many  types  of  molecular  excitations  can  be  de¬ 
scribed  by  relatively  few  equations.  For  optically  allowed  excitations, 
the  integrated  cross  section  is  given  by 
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II 


4?ra^H 


X  WC„E) 


(44) 


which  is  identical  in  form  to  Eq.  (31)  for  atomic  processes.  For  ex¬ 
citation  to  a  stable  molecular  state  n,  is  related  to  the  optical 
oscillator  strength  fn  by  =  fn(R/En)  where  En  is  the  excitation 
energy.  In  the  case  of  dissociative  excitation  to  the  atomic  state  n, 
is  related  to  the  optical  oscillator  strengths  of  all  dipole  transitions 
to  molecular  states  which  lead  to  formation  of  the  particular  state  n  and 
is  given  by 


n 


(45) 


where  E'  is  the  excitation  energy  transferred  to  the  molecule,  En  is 
the  threshold  energy  for  the  process  under,  consideration,  r}n(E')  is  the 
probability  that  the  intermediate  molecular  states  will  dissociate  into 
the  particular  atomic  state  n,  and  df/dE"  is  the  differential  optical 
oscillator  strength.  The  remaining  coefficient  Cn  in  the  cross-section 
formula  is  obviously  defined  by 


M„Fn  C 
q  n 


?n  C  E  /dE  ' 


(46) 


For  ionization,  the  factor  rjn(E')  need  only  be  redefined  as  the  probabil¬ 
ity  that  the  molecule  ionizes  upon  receiving  an  energy  transfer  E". 


For  optically  forbidden  excitations,  the  integrated  cross  section  is 
given  by 


- —  b 

t  n  (47) 


where  bn  is  a  constant  related  by  Eq.  (33)  to  the  generalized  oscillator 
strength  defined  in  terms  of  the  molecular  form  factor.  In  the  case  of 
dissociative  excitation  to  the  atomic  state  n,  bn  is  a  constant  depending 
on  the  sum  of  all  optically  forbidden  processes  leading  to  formation  of 
atoms  in  state  n. 

Whenever  the  cross  sections  are  plotted  in  the  form  oE/(47raoR) 
versus  in  E(Bethe  plot),  the  optically  allowed,  optically  forbidden,  and 
spin  forbidden  excitations  have  curves  with  positive,  zero,  or  negative 
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2 

slope,  respectively.  The  slopes  of  these  lines  are  related  to  Mn  and 
the  intercepts  to  Cn.  It  is  generally  found  that  for  optically  allowed 
processes  Cn  has  a  value  near  one  and  for  optically  forbidden  processes 
Cn  »  1  (Ref.  4). 

In  summary,  to  enable  the  knowledgeable  design  of  electron  beam 
systems  for  investigations  of  flow  fields  of  such  binary  mixtures  as 
N2/He  and  N2/NO,  it  is  desired  to  experimentally  determine  the  func¬ 
tional  dependence  on  energy  (see  Eqs.  (44)  and  (47))  of  various  dissoci¬ 
ative  excitation  processes  of  electron-N2  and  -NO  collisions.  Further, 
the  variation  of  the  electron-beam-induced  fluorescence  intensity  with 
both  beam  current  and  molecular  specie  density  is  to  be  investigated 
under  such  experimental  conditions  as  to  preclude  collisional  quenching 
effects.  The  experimental  apparatus  used  for  these  purposes  is  de¬ 
scribed  in  the  following  section. 


2.0  EXPERIMENTAL  APPARATUS 


2.1  VACUUM  CHAMBER  AND  GAS  SUPPLY 

The  vacuum  chamber  was  a  stainless  steel  cylindrical  vessel  of 
1  26.  5  cm  inside  diameter  and  45  cm  length.  Interior  to  the  vacuum 
chamber  was  installed  a  cylindrical  copper  cooling  liner  which  was 
aligned  both  axially  and  radially  with  the  vacuum  chamber.  The  cooling 
liner  was  of  22  cm  inside  diameter,  extended  the  length  of  the  vacuum 
chamber,  and  was  held  in  place  by  several  set  screws.  The  cooling 
liner  had  three  2.  5-cm-diam  access  ports,  one  of  which  was  used  for 
the  electron-beam  injection  orifice.  The  remaining  two  ports  were 
available  for  optical  detection  of  the  electron-beam  fluorescence  and 
were  aligned  with  the  quartz  optical  window  of  the  vacuum  chamber. 

The  vacuum  chamber  was  essentially  that  of  Ref.  5. 

The  copper  liner  temperature  was  determined  by  a  copper- 
constantan  thermocouple  attached  to  the  liner,  and  the  observed  vari¬ 
ation  of  liner  temperature  was  less  than  0.  5  K. 

Chamber  pumping  was  provided  by  a  10.  2-cm-diam  oil  diffusion 
pump  with  a  water-cooled  baffle,  providing  chamber  pressures  of  less 
than  10 “6  torr.  The  baffle  cooling  source  was  later  changed  from  water 
to  a  Freon®  refrigerator  to  further  reduce  diffusion  pump  backstr earning. 
The  chamber  pressure  was  controlled  by  varying  the  mass  flow  rate 
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through  the  system  using  both  a  needle  valve  and  a  Granville- Phillips 
leak  valve  at  the  inlet  of  the  gas  to  the  chamber  and  also  by  adjusting 
the  gate  valve  opening  to  the  10.  2 -cm  diffusion  pump  at  the  chamber 
exit. 


An  NRC  Model  820  Alphatron®  gage  monitored  the  pressure  at  a 
point  between  the  copper  liner  and  chamber  wall;  a  calibrated  MKS 
Baratron  was  used  to  determine  the  scattering  cell  pressure  for  the 
relative  intensity  measurements. 

The  mechanical  pump  exhaust  was  vented  to  the  atmosphere  through 
a  1.  2  7-cm-diam  copper  tube.  To  avoid  back  pressure  on  the  mechanical 
pump  and  to  dilute  the  concentration  of  NO  in  the  vent  tube,  an  injector 
nozzle  was  attached  to  the  pump  vent  connection,  and  high-pressure  dry 
nitrogen  was  used  to  purge  the  venting  tube. 


2.2  ELECTRON  GUN  SYSTEM 

The  electron  gun,  which  was  the  same  as  used  for  the  measurements 
reported  in  Ref.  5,  was  a  television  type,  oxide-coated  cathode  model 
and  typically  provided  a  steady  1.5-  to  3. 0-mA  current  for  10  to  30  hr. 

A  Universal  Voltronics  high-voltage  power  supply  with  a  0.01-percent 
ripple  filter  provided  a  negative  potential  of  4  to  20  kV  to  the  electron 
gun  cathode.  The  high  voltage  power  supply,  which  was  calibrated 
using  a  secondary  standard  Weston  ammeter  and  calibrated  resistors, 
had  an  error  of  less  than  5  percent  between  4  and  7  kV  and  less  than 
2  percent  above  7  kV  as  shown  in  Fig.  4.  The  filament  and  accelerator 
power  supply  shown  in  Fig.  5  was  designed  to  replace  batteries  used 
previously  and  to  allow  immersion  in  oil,  which  permitted  operation  up 
to  50  kV  without  arcing  and  corona  effects. 

Magnetic  coils  situated  near  the  electron  gun  cathode  provided  for 
electron  beam  focus  at  and  deflection  about  the  1.  3-mm-diam  orifice, 
which  was  situated  approximately  56  cm  from  the  gun  cathode  and 
located  just  interior  to  the  vacuum  chamber  as  shown  in  Fig.  6.  A 
10.  2 -cm  oil  diffusion  pump  with  a  water-cooled,  and  later  Freon- 
refrigerator-cooled,  baffle  maintained  the  differentially  pumped  drift 
tube  section  at  a  pressure  of  less  than  10~5  torr  when  the  chamber  pres¬ 
sure  was  2  x  10“3  torr. 
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Power  Supply,  kv 

Figure  4.  Calibration  of  Universal  Voltronics  high  voltage 
power  supply. 


Figure  5.  Electron  gun  filament  and  accelerator  power  supply. 
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The  electron  beam  Faraday  cup  collector  consisted  of  a  water- 
cooled,  4.  1-cm-diam  right-angle  elbow  copper  cylinder  and  two  screen 
grids  with  1.  9-cm-diam  entrance  holes,  the  second  of  which  was  biased 
at  -30  v  for  the  suppression  of  secondary  electrons.  The  initial  collect¬ 
ing  surface  was  inclined  at  an  angle  of  45  deg  to  the  beam  direction  and, 
therefore,  extended  from  2.  5  to  4.  1  cm  below  the  top  grid.  Measuring 
from  the  center  of  the  inclined  surface,  the  cylinder  extended  12.  7  cm 
to  collect  scattered  electrons.  This  particular  configuration  of  the 
Faraday  cup  was  necessary  to  provide,  within  the  space  limitations  of 
the  chamber,  an  adequate  length-to-diameter  ratio  of  the  collector  to 
minimize  reflective  loss  of  the  electron  beam  on  inpact  with  the  collector. 
The  interior  surfaces  of  the  Faraday  cup  were  coated  with  Aquadag® ,  a 
colloidal  graphite  suspension.  The  beam  current  collected  by  the 
Faraday  cup  was  determined  by  measuring  the  potential  drop  across  a 
10-kf2  resistor  with  an  NLS  Model  481  digital  voltmeter,  which  was  later 
replaced  with  a  Heathkit  Model  IM-102  multimeter. 


2.3  OPTICS  AND  DETECTORS 

The  optical  arrangement  employed  for  these  measurements  is  shown 
in  Fig.  6.  The  10-cm-diam  glass  lens  had  a  focal  length  of  25.  54  cm, 
and  the  overall  optical  magnification  of  the  system  was  1.0.  AO.  5-m 
Jarrell  Ash  scanning  spectrometer  was  employed  with  an  1180-groove/mm 
grating  blazed  at  5000  A,  and  the  reciprocal  linear  dispersion  was 
16  A/mm.  The  spectrometer  was  mounted  on  its  side  with  the  entrance 
slit  of  the  spectrometer  imaged  onto  the  chamber  centerline  such  that  the 
slit  length  was  perpendicular  to  the  electron  beam  direction. 

An  EMI  6256S  photomultiplier  tube  of  S-ll  spectral  response  was 
used,  and  the  operating  voltage  of  1500  v  was  supplied  by  a  Fluke  Model 
405B,  high-voltage  power  supply.  The  photomultiplier  tube  tempera¬ 
ture  was  maintained  near  -32°C,  using  gaseous  N2  cooling,  for  the  pur¬ 
poses  of  reducing  the  dark  count  and  increasing  the  overall  signal -to- 
noise  ratio  of  the  tube  (Refs.  6  and  7).  The  dark  count  was  normally  one 
count  per  second  at  this  operating  temperature. 

The  relative  spectral  sensitivity  of  the  system  was  determined  by 
placing  a  GE  30A-T24-17  tungsten  strip  lamp  at  the  observation  volume 
and  recording  the  system  response  as  a  function  of  wavelength.  The 
spectral  radiance  N\(X,  Tg)  of  the  lamp  operating  at  a  true  lamp  tem¬ 
perature  T  with  emissivity  e(\,  T)  and  lamp  transmissivity  t  is  given  by 
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\^(A.'I  |j)  =  r  *(A,T)  \^(A,T).  watts  '  cm^-sr-cm 

or  in  terms  of  the  photon  rate 

VA.TV  =  N\(A.Tu)'(hA'c) 


(48) 


where  A.  is  the  wavelength  of  the  radiation  and  Tg  is  the  brightness 
temperature  of  the  lamp.  It  is  seen  from  Eq.  (48)  that  Tg  is  the  tem¬ 
perature  of  a  blackbody  which  emits  the  same  radiant  intensity  at  X  as 
does  the  tungsten  surface  at  temperature  T. 

Using  the  Planck  radiation  law  the  spectral  radiance  is  given  by 


Xa(A.Tb)  =  C jA‘S(exp  (Cj'ATy)  -  l]1.  watts  cm2-sr-cm 


(49) 


By  assuming  large  values  of  C2/A.T23,  one  can  write  Eq.  (49)  in  the 
Wien  law  approximation  as 


VA.TV  =  C]A'5  exp  [-C2/ATg]  (50) 

Consequently,  Eq.  (48)  becomes 

C|A"S  exp  (-C2  ’ATb)  =  r  e(A.T)  ■  C  |A'5[exp  (C2'  AT)  -  l]1  (51) 

By  using  the  Wien  law  approximation,  Eq.  (51)  can  be  written  as 

l  'T  =  1  Tjj  +  (A'  C2)  fnMA.T)]  (52) 

The  emissivity  data  of  DeVos  (Ref.  8)  show  that,  for  X  =  6500  A,  the 
pyrometer  wavelength  e(650  nm,  T)  can  be  written  as 


e(650  nm,T)  =  <  (650  nm.2400)  +  [(T  -  2200)  "200][f  (650  mri,22nn)  -  e (650  nm. 2400)1  (53) 

for  2200  <  T  <  2400  K.  Since  t  for  the  fused  silica  strip  lamp  window 
is  known  to  be  0.  93  (Ref.  9),  Eqs.  (52)  and  (53)  can  be  solved  itera¬ 
tively  for  the  true  lamp  temperature  T  once  the  brightness  temperature 
is  determined  using  a  pyrometer.  Using  the  true  temperature  T  pro¬ 
duces  the  spectral  radiance  of  the  lamp  of 
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Nx  =  r  e(A,T)  •  CjA-5[exp(C2/AT)  -  l]’1 


(54) 


The  constants  Ci  and  C2  are 


Cj  =  1.1909  x  10-12,  vat t- cm2/ sr 


C2  =  1.4380  ,  cm-K 


The  relative  spectral  sensitivity  S(X)  of  the  system  is  defined  as 
the  ratio  of  the  system  response  normalized  by  the  lamp  radiance 

Nx  and  the  maximum  value  of  for  the  spectral  range  studied; 


1.  e. 


(55) 


The  relative  spectral  sensitivity  S(X)  was  determined  by  performing  rela¬ 
tive  intensity  measurements  of  the  lamp  output  at  50-A  intervals  over  the 
spectral  range  of  3500  to  6550  A  using  200-Mm  spectrometer  slit  widths 
and  two  neutral  density  filters  of  known  transmission  factors.  The  strip 
lamp  was  powered  by  an  EG  and  G  Model  590-11  power  supply  operated 
at  30  amp  and  6  v.  This  supply  provides  a  chopper-stabilized,  1-kHz 
square-wave  output  which  is  feedback  regulated  to  within  0.25  percent 
rms  of  the  selected  value  which  provides  control  of  the  lamp  output  to 
within  1  percent.  A  calibrated  L  and  N  optical  pyrometer  was  used  to 
determine  the  brightness  temperature  Tg  of  the  lamp  which  was  2141  ±6  K 
for  these  calibrations.  An  iterative  computer  calculation  was  used  for 
Eqs.  (52)  to  (54)  to  determine  S(A)  given  by  Eq.  (55),  and  S(A)  for  the 
system  as  shown  in  Fig.  7.  It  is  estimated  that  the  relative  S(A)  values 
are  determined  to  within  an  inaccuracy  of  1  to  2  percent. 


2.4  PHOTON  COUNTING  SYSTEM 

The  photon  counting  system,  which  was  the  same  as  that  reported 
in  Ref.  5,  consisted  of  the  following  ORTEC  apparatus:  Model  454 
Timing  Filter  Amplifier,  Model  436  100  MHz  Discriminator,  and 
Model  715  Dual  Counter/ Timer.  Figure  8  shows  a  block  diagram  of  the 
system. 
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Figure  7.  Relative  spectral  response  of  the  system. 
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Figure  8.  Block  diagram  of  ORTEC  photon  counting  system. 
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The  715  counter /timer  includes  dual  counters,  one  of  which,  in 
conjunction  with  a  precise  internal  oscillator,  may  be  preset  for  an 
accurate  counting  interval  so  that  the  system  can  be  used  as  a  discrete 
rate  meter.  The  output  is  available  as  either  positive  or  negative 
pulses.  The  negative  output  is  passively  shaped  to  an  extent  deter¬ 
mined  by  the  time  constant  selection,  but  it  is  required  only  for  ex¬ 
tremely  high  rate  applications.  The  positive  pulse  is  a  standardized 
pulse  of  +5  v  and  0.  5  jusec  width.  The  715  counter  has  a  frequency 
response  of  20  MHz  and  an  approximate  dead -time  of  50  nsec.  The 
counter  /timer  will  accept  either  discriminator  output,  and  for  this  re¬ 
port,  the  negative  output  was  used  for  the  measurements  that  were 
made.  The  discriminator  was  operated  at  0.  6  v  for  the  purpose  of  ex¬ 
cluding  dark  counts  and  obtaining  the  highest  signal-to-noise  ratio  of 
the  system  (Ref.  10). 


3.0  RESULTS  AND  ANALYSIS 


Preliminary  to  the  cross-section  measurements,  atomic  line  spec¬ 
tral  profile  data  wereo acquired  using  100-f^m  slit  widths  giving  a  spec¬ 
tral  bandpass  of  1.6  A  in  order  to  isolate  spectral  impurities  and  deter¬ 
mine  appropriate  slit  widths.  Since  the  slit  length  was  perpendicular  to 
the  beam,  the  photon  count  rate  was  proportional  to  the  slit  width,  and 
increasing  the  slit  width  pays  handsome  dividends.  The  spectral  line 
profiles  were  determined  by  measuring  the  radiative  photon  rate  of  each 
transition  as  a  function  of  wavelength.  Figures  9  through  12  show  the 
spectral  feature^  of  NIKN2;  5016.4  A),  NII(NO;  5016.4  A), 

NII(N2;  3995.0  A),  and  NII(NO;  3995.0  A),  respectively.  The  data  were 
acquired  at  a  beam  energy  of  10  keV,  except  for  NII(N2;  3995.  0  A)  which 
was  at  7  keV,  and  the  pressures  and  currents  were:  8.35  mtorr,  860  juA; 
5.  3  mtorr,  745  pA;  2.  8  mtorr,  220  M A;  and  5.  5  mtorr,  725  iu A,  respec¬ 
tively. 

The  NII(N2;o3995.  0  A)  spectral  impurity  is  due  to  the  N2(2+)(l,4) 
band  at  3998.  4  A  (Ref.  11).  By  using  the  peak  cross  sections  for  the 
N2(2+)(l,4)  band  (Ref.  12),  the  secondary  electron  contribution  to  the 
NIKN2;  3995.  0  A)  current  was  found  to  be  less  than  0.  05  percent  over 
the  4  to  21  keV  electron  beam  energy  range. 

By  using  the  NII(NO;  3995.  0  A)  line,  which  appears  to  be  free  of  all 
spectral  impurities,  the  full  width  at  half  maximum  is  found  to  be  1.6  A 
in  agreement  with  the  assumed  reciprocal  linear  dispersion  of  16  A/mm 
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and  slit  width  of  100  jum.  A  slit  width  of  200  jum  was  used  for  the  cross- 
section  measurements  since  the  observed  intensity  was  proportional  to 
the  slit  width  and  200  jum  was  still  sufficiently  small  to  safely  resolve 
the  5016. 4 -A  lines. 

The  current  and  pressure  dependence  was  determined  by  measur¬ 
ing  the  radiative  photon  rate  of  each  transition  as  a  function  of  the  beam 
current  and  chamber  pressure,  respectively.  Figures  13  to  16  show  the 
features  of  the  current  and  pressurg  dependence  of  NII(N2;  5016.4  A), 
NH(NO;  5016.4  A),  NII(N2;  3995.0  A),  and  NII(NO;  3995.0  A),  respec¬ 
tively.  The  data  were  acquired  at  a  beam  energy  of  10  keV  and  a  slit 
width  of  200  jum  except  for  the  pressure  dependence  of  NII(N2;  3995.0  A), 
which  was  obtained  at  5  keV,  and  the  NII(NO;  3995.  0  A)  current  depend¬ 
ence,  which  was  obtained  at  15  keV.  In  each  case,  the  pressure  was 
maintained  at  a  constant  value  for  the  current  dependence  measurements, 
and  the  beam  current  was  maintained  at  a  constant  value  for  the  pressure 
dependence  measurements. 
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Figure  9.  NH(N2;  5016.4  A)  line  profile  (slit  width  =100  jum). 
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Figure  10.  NII(NO;  5016.4  A)  line  profile  (slit  width  =  100  /um). 
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Figure  13.  Current  and  pressure  dependence  of  the  molecular 
nitrogen  Nll(3  F°  -*■  3D)  line  at  5016.4  A. 
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Figure  14.  Current  and  pressure  dependence  of  the  nitric  oxide 
NIK3  F°  3D)  line  at  5016.4  A. 
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Figure  15.  Current  and  pressure  dependence  of  the  molecular 
nitrogen  NIM1  D  1  P°)  line  at  3995.0  A. 
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Figure  16.  Current  and  pressure  dependence  of  the  nitric 
oxide  Nil t1  D  -»■  1  P°)  line  at  3995.0  A. 
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The  current  and  pressure  dependences  were  obtained  so  that  the 
cross-section  measurements  could  be  made  at  currents  and  pressures 
minimizing,  or  free  of,  secondary  electron  and  collisional  quenching 
effects  and  sufficiently  linear  to  permit  normalizing  to  the  current  and 
pressure.  Figures  13  and  14  show  the  variations  of  the  intensities  for 
the  two  5016.  4-A  transitions  with  both  current  and  pressure  to  be  linear 
over  the  observed  range  and  to  pass  through  the  origin,  thus  allowing 
the  desired  normalization.  Similarly,  Figs.  15  and  16  show  a  linear 
current  dependence  for  both  3995.  0-A  transitions  over  the  observed 
range  of  beam  current  and  the  intensities  extrapolate  to  zero  for  zero 
beam  current.  The  pressure  dependence  of  the  NIKN2;  3995.  0  A)  line 
is  linear  over  the  observed  pressure  range  but  does  not  pass  through 
the  origin,  while  the  NII(NO;  2995.  0  A)  line  exhibits  definite  non¬ 
linearity  at  the  lower  pressures.  In  both  cases,  however,  the  linear 
region  at  high  pressures  can  be  extrapolated  to  an  abscissa  intercept 
of  about  0.  7  mtorr. 

The  energy  dependence  of  the  optical  excitation  cross  sections  was 
determined  by  measuring  the  radiative  photon  rate  of  each  transition  as 
a  function  of  the  electron  beam  energy.  By  using  a  spectrometer  slit 
width  of  200  Aim,  the  total  photon  count  summation  at  each  wavelength 
of  interest  was  accumulated  over  the  time  interval,  less  the  photomulti¬ 
plier  tube  dark  count,  was  divided  by  the  beam  current,  chamber  pres¬ 
sure,  and  count  interval  time.  The  operating  conditions  were  chosen 
such  that,  the  linearity  of  intensity  with  beam  current  and  chamber  pres¬ 
sure  was  ensured  (Table  1). 

Table  1.  Experimental  Conditions  during  Data  Acquisition 


Specie 

Energy  Dependence 

Ratios  at  10  keV 

Beam 

Current  (pA) 

Pressure 

(mtorr) 

Beam 

Current  (p  A) 

Pressure 

(mtorr) 

N2+(H  (0,0) 

<825 

<3.0 

425 

3.0 

Nil  <N2;  5016. 4  A) 

<800 

10.2 

425 

8.7 

Nil  %  3995.0  A) 

<240 

2.7 

285 

3.0 

Nil  (NO;  5016.4  A) 

<995 

6.0 

355 

3.8 

NIMN0;  3995.0  A) 

<995 

6.6 

355 

3.4 
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For  the  purpose  of  determining  absolute  values  of  the  excitation 
cross  section  without  the  requirement  of  making  absolute  intensity  mea¬ 
surements,  it  was  decided  that  measurements  of  the  strong,  optically 
allowed  N^CO.O)  band  of  the  N2  First  Negative  System,  Njjfl-),  would  be 
performed  and  the  results  used  for  normalization  of  the  Nil  results. 

This  particular  radiative  system  was  selected  because  it  has  been  veri¬ 
fied  that  it  obeys  the  prediction  of  the  Bethe-Born  relation  for  its  energy 
dependence  and  absolute  excitation  cross  sections  are  available  for  elec¬ 
tron  impact  energies  up  to  6  keV  (Ref.  13).  Additionally,  adequate  theo¬ 
retical  and  experimental  knowledge  exist  for  the  collision-radiation  pro¬ 
cess 

NjX'SJv  =  0)  +  e-  ->  NljB2£+(v  =  0)  +  2c-  (56) 

L  NtX22+(v  =  0)  +  hi. 


to  enable  accurate  calculations  of  the  intensity  and  profile  of  the  (0,0) 
band  structure  over  a  very  wide  energy  range  (Ref.  14).  The  energy- 
dependence  of  the  partial  N^U-HO,  0)  band  intensity  was  obtained  by 
setting  the  spectrometer  bandpass  of  3.2  A  on  the  P  branch  peak  and 
measuring  the  photon  rate  as  a  function  of  the  electron  energy.  The 
data  of  Ref.  13  were  used  to  obtain  an  extrapolated  value  of  9.  93  x  10" 19 
cm^  for  the  N2(l-)(0,0)  band  absolute  emission  cross  section  at  the  elec¬ 
tron  beam  energy  of  10  keV,  and  the  experimental  results  of  the  present 
study  were  normalized  to  this  point.  The  results  are  shown  in  Fig.  17, 


Figure  17.  Bethe  plot  of  the  N£(1-)  (0,0)  band  optical 
excitation  cross  section. 
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o 

and  it  is  seen  that  the  measured  values  of  aE/4?raoR  are  linear  with 
j£n  E,  verifying  that  the  excitation  of  the  N^U-HO,  0)  band  is  optically 
allowed  over  the  energy  range  studied.  These  results  are  shown  as 
well  in  Table  2.  The  slope  obtained  for  the  present  results  is 
0.  353  ±  0.  009  with  an  intercept  of  -1.  18  ±  0.  08  in  agreement  with  a 

Table  2.  Energy  Dependence  of  the  N£(1-)  (0,0)  Band 
Emission  Cross  Section 


E(keV) 

oE 

o(10"^  cm^) 

4Tra0^R 

4.0 

1.75 

20.9 

6.0 

1.88 

15.0 

10.0 

2.08 

9.96 

15.0 

2.21 

7.05 

20.0 

2.33 

5.58 

slope  of  0.  35  given  in  Ref.  13.  The  standard  deviation  lis,ted  is  the 
one -sigma  value.  Figures  18  and  19  show  the  current  and  pressure 
dependence  measurements  made  at  a  beam  energy  of  10  keV  and  a  slit 
width  of  200  dm  and  the  pressure  dependence  at  4  keV,  respectively. 
In  each  case  the  observed  radiative  photon  rate  is  linear  and  passes 
through  the  origin  permitting  the  desired  normalization  to  current  .and 
pressure. 


Figure  18.  Current  and  pressure  dependence  of  the  NJO-)  (0,0)  band. 
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Pressure  (mtorr) 

Figure  19.  Pressure  dependence  of  die  NJ(I-)  (0,0)  band  at  a 
beam  energy  of  4  keV. 

To  use  the  NjjU-XO,  0)  band  intensity  for  placing  the  Nil  data  on  an 
absolute  scale,  the  entire  band  must  be  measured.  The  maximum  slit 
width  of  the  0.  5-m  spectrometer  used  was  400  Mm  giving  a  bandpass  of 
approximately  6.  4  A,  which  is  less  than  approximately  1/4  that  of  the 
bandwidth.  Consequently,  it  was  necessary  to  adjust  the  measured 
partial  band  intensity  to  obtain  the  required  total  band  intensity  value. 

The  equations  required  for  the  theoretical  calculation  of  the  spectral 
intensity  of  a  vibrational  band  in  the  N^U-)  system  are  given  by  Williams 
(Ref.  14).  These  equations  enable  the  calculation  of  the  theoretical  rela¬ 
tive  intensity  of  each  rotational  line  of  the  band  as  a  function  of  the  vi¬ 
brational  and  rotational  temperatures,  assuming  dipole  rotational  excita¬ 
tions  and  including  nuclear  spin  considerations.  The  set  of  emission 
band  strengths,  excitation  Franck-Condon  factors,  and  other  molecular 
parameters  used  in  a  computer  calculation  of  the  1^(1 -)  system  rotational 
line  intensities  are  listed  or  referenced  in  Ref.  5.  The  rotational  line 
wavelength  calculation,  contained  in  the  computer  calculation,  included 
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corrections  for  anharmonicity  of  the  vibration  and  the  influence  of  cen¬ 
trifugal  stretching.  The  computer  calculation  then  convolved  the  theo¬ 
retical  rotational  line  intensities  with  the  spectrometer  response  func¬ 
tion,  using  the  reciprocal  linear  dispersion  and  spectrometer  slit 
widths  as  input  data,  to  calculate  the  theoretical  spectrometer  response 
from  a  given  vibrational  band  of  the  N^U-)  system  excited  by  an  elec¬ 
tron  beam  in  0.25-A  increments.  The  resulting  plot  of  the  Nr^l-MO.O) 
band  at  280  K  and  200  Mm  slit  width  is  shown  in  Fig.  20  and  may  be  com¬ 
pared  with  the  experimentally  obtained  profile  in  Fig.  21.  The  theoreti¬ 
cal  ratio  of  the  P  branch  peak  at  a  slit  width  of  200  Mm  to  the  entire  band 
was  calculated  to  be  0.423.  Variations  of  0.  25  Ain  the  spectrometer 
position  near  the  P  branch  peak  were  shown  by  the  calculation  to  change 
the  partial  band  intensity  by  approximately  0.  1  percent  of  the  total  band 
intensity. 

The  intensities  of  the  N2<l-K0,0)  band  and  the  Nil  lines  were  obtained 
using  a  slit  width  of  200  Mm  at  an  electron  beam  energy  of  10  keV.  The 
measured  radiative  photon  count  rates  were  corrected  for  photomultiplier 
tube  dark  count,  and  then  divided  by  the  beam  current  and  chamber  pres¬ 
sure.  The  results  were  corrected  for  background  emission  (5016.  4-A 
lines)  and  then  divided  by  the  relative  spectral  sensitivity  factor  S(A)  to 
obtain  relative  intensityo measurements.  The  background  correction 
factors  for  the  5016.  4-A  Nil  lines  were  obtained  by  subtracting  the  ratio 
of  the  background  intensity  and  the  total  line  intensity  from  unity.  '  The 
relative  intensity  measurements  for  each  spectral  line  were  multiplied 
by  the  calculated  partial  to  entire  NjjU-MO,  0)  band  ratio  and  the 
Nj-jU-XO,  0)  band  absolute  emission  cross  section  to  obtain  the  absolute 
optical  excitation  cross  sections  and  subsequently  the  apparent  cross 
sections.  Table  3  shows  the  results  of  the  measurements  at  the  electron 
impact  energy  of  10  keV. 

For  the  purpose  of  analyzing  the  numerous  energy  scans  made  for 
each  transition,  a  computer  program  was  developed  and  is  shown  in 
Table  4.  The  photon  count  rate,  photomultiplier  tube  dark  count,  beam 
current,  counting  times  for  the  intensity  and  dark  count,  and  chamber 
pressure  at  each  electron  beam  energy  were  used  as  input  data.  Also 
input  was  the  normalized  value  of  the  ordinate  at  an  electron  beam 
energy  of  10  keV,  which  is  simply  the  absolute  apparent  cross  section 
for  the  plot  using  the  reciprocal  beam  energy  as  the  abscissa  ancfiiP 
equal  to  crE/4jraQR  for  the  Bethe  plot  using  i n  E  as  the  abscissa. 
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Wavelength,  angstroms 


Figure  20.  Computer  calculated  N£(1-}  (0,0)  band  profile. 


Figure  21.  N^M-)  (0.0)  band  profile  (slit  width  =  200  nm). 
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Table  3.  Nil  Apparent  Cross  Sections  at  10  keV. 


Specie 

x,  A 

Array 

Multiple! 

Background 

Correction 

Factor 

o(X) 104 
0(0,0) 

Optical* 

o(X) 

(10“22  cm2) 

Pui 

Apparent 

o(X) 

(10“21  cm2) 

OE103 
47ra02  R 

D 

5016.4 

2p3d-**2p3p 

3f°— 3d 

a  735 

2.25 

2.23 

a  155 

1.44 

3.01 

3995.0 

2p3p-*-2p3s 

lp-^lpo 

1 

2L6 

21.4 

i 

2.14 

447 

NO 

5016.4 

2p3d— 2p3p 

3(io_^3D 

0.753 

1.94 

1.93 

a  155 

L25 

2.61 

3995.0 

2p3p-*-2p3s 

l[)_*.lpO 

1 

17.7 

17.6 

i 

1. 76 

3.68 

•Normalized  to  the  N2+(l-)  (0, 0)  Band  Emission  Cross  Section  at  10  keV,  9. 93  x  10"19  cm2 
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0001 

0002 


0003 

0004 

0005 

0006 

0007 

0006 

000.9 

0010 

0011 

0012 

0013 

0014 

0015 

0016 

0017 

0016 

0019 

0020 

0021 

0022 

0023 

0024 

0025 

0026 


0027 

0026 


0029 


0030 


Table  4.  Program  SUMDATA 


C  HRuGWAM  SUMO A T A 

DIMENSION  bE (50) »CH (30) (UK6 (3U) »BC(30> ,TCH (30 ) . TbKG <30 ) * 

1  P ( 30  )  . lHUW (400) * I COL (400 ) i 1 COUNT (20) 

UUUBLt  PRECISION  SCR.USCH.S130) »D(3U) .AVSU.wT (30) .FU6t(5Q) t 

1  5W  I  » b*  I  A  *  SW  f  XX  t  S*l  T  Y  »  SW  TYY »  SwTXY  <  SWMl  ,SWWTX»SWWTXX.UtN(2U> • 

2  SLP (20) .CEP (20) .SY(20> ,D5LP(20) .DLtP(20> .REFE0.RNORM120) t 

3  UElHt (50 ) »DBET  HE ( 50) .BEX(20.2U) .BETHX (2U .20 ) tUBETMX (20 »20) • 

4  ABt(4U0) * XbETM ( 400 ) .XUUETH (400) ♦  A »CXBE  (4|)Q ) fSUM.SSUM, 

5  ABETHE (50) « AUbETH (50) »l)EM , SlM .CEM i SM.OSLM.OCEM 

400  FUHMAl <10X,F10.U*6(5X,1PE12.S)/I 

401  FORMAT I5X.4(5XtlPE12*5)/> 

402  FUHMAl  (6A,2(1HE12.5) ) 

501  FORMAT  (13) 

502  FURMAI  (6X, 1PE12.5. 17) 

503  FUHMAl  (6A«F10.3»2F10.2»2F10.0»F10.3) 

504  FORMA!  (20X»5(5X* 1PE12.5) ./) 

507  FUHMAl  (315) 

550  FUHMAl  (//.25X*9H  SLOPE  s  , lPtl2.5»3X* 16H  WITH  STD.  DEV,  ■  > 

1  1PE12.5./) 

551  FUHMAK//.25X.13H  INTERCEPT  *  .IPE12.S.  .  3*.  _ _ 

1  lOH  Hi f M  510.  DEV.  «  • lPt!2.5) 

603  FUHMAl (///.35X.16H**  BE THE  PLOT  •*•///) 

604  FUHMAl (///.35X.26H**  RECIPROCAL  tNERGY  PLOT  •«»///) 

605  FUHMAl (// »b5X. 19HUATA  NORMALIZED  AT  .F7.0.3H  EV«31H  USING  AN  OR 
1D1NATE  EUUAL  TO  . 1PE12.5./// > 

900  FORMAT  (20A4)  ........  -  — 

901  FUHMAl  I1H1.25A.2UA4.//) 

902  FUHMAl (//.25X.20A4.//) 

903  FUHMAl  (//.25X.20A4.//1 

1000  FUHMAl  <//ilOX.'  ENtRGY  A66C1SSA  ORDINATE 

1  STD.  OEV.  HEIGHT  BETHE  STD.  DEV.'t/) 

.  1010  FUHMAl  (//.10X»»  SLOPES  StU.  DEV..  INTERCEPTS...  . 

1  STU.  UbV.  «./> 

1020  FORMAT (//»2BX»6HEN£RGY»1QA*8HAbSCISSA.9X»6H0RDINATE»9X» 

1  9HSTD.  DEV.  .9X.6HVIE1GHT.//) 

READ  900f A1 .A2. A3.A4.A5. A6.A7.A6.A9.A10tAll.A12.A13.A14.A15> 
1A16.A17.A16.A19.A2U 

_  READ  900.61 .U2.63.U4.B5.B6.U7. 66. B9.B10. 611. 612. H13.bl4.B15t, _ 

lB16.Bi7.B16.B19.H2U 

READ  900iCl.C2iC3.C4.C5.C6.C7.C6.C9.C10.Cll.C12.C13.C14.C15. 

1  C16.L17.Cl6fC19.C20 

READ  402.SIGMA0.E0 

C  S1GMA0  15  THE  AB5ULUTE  CROSS  SECTION  AT  ENERGY  EO  FOR  A 

C  RECIPRULAL  ENEHGY  PLOT.  .  .  _  -  ...... 

C  -  SiGMAO  15  THE  OHUINATE  OF  A  BETHE  PLUl  AT  ENERGY  EO  WHEN  JG»1 
JG  =  1 

READ  507.  1CMAX. IRMAX.NVAL 

C  IHMAX  is  The  number  OF  ENERGY  POINTS  In  the  largest  daia  set 

c  ilmax  is  the  number  uf  sets  of  data 

C  NVAL  is  the  PRODUCT  ICMAX*IRMAa  AND  DEFINES  THE  NUMBER  OF 

C  ELEMENTS  IN  THE  ARRAY 

5  HEAD  501. NR 

C  NA  IS  THE  NUMBER  UF  SETS  OF  DATA.  OR  THE  NUMBER  OF  SLUPES. 

L  hhich  will  be  calculated 

PRINT  90l.Al.A2.A3.A4.AStAb.A7.A6.A9.A10.All.A12.A13.A14.AlS. 
1A16.A17.A16.A19.A20  .... 
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Table  4.  Continued 


0031 

0032 

0033 

0034 

0035 

0036 

0037 

0030 

0039 


0040 

0041 


0042 

0043 


0044 

0045 

0046 

0047 

0040 

0049 

0050 


0051 

0Ub2 


00b3 

0054 

0055 

0056 

0057 

0050 

0059 

0060 


OObl 

0062 

00o3 

0064 

0065 


PRINT  902, Bl, 02. 03, B4, b5, 06. b7,bb, 09. BIO, 011,012.013,014.015. 
Iul6*0l7.0iefbl9f020  .... 

POINT  903»Cl,C2»C3,C4,C5,C6»CT»CB»C9»C10«Cll,C12iC13»C14»Cl5» 

1  C lb, Cl 7 .Clb.tl9.C20 

1FUG.EU.2)  GO  TO  4 

3  PRINT  603 
(,0  TO  9 

4  PRINT  604 

9  CONTINUE 

DU  65  Kal»NK 
10  HEAD  50 1 • NO 

C  NJ  IS  THE  NUMBER  UF  ENEHGIES  mhERE  DA  I A  NAS  TAKEN  TO  OBTAIN  A 
C  SET  SUFFICIENT  TU  DETERMINE  THE  5LUPE 
DU  30  Jal,NJ 
15  RLaO  502.0E ( J1 »Ni 

C  bt  IS  THE  OEAM  ENEHGY  IN  tLECIRUN  VOLTS 

C  hi  IS  THt  NUMBER  OF  VALUES  TAKEN  AT  EACH  ENERGY 

DU  20  Iai *Ni 

20  READ  503*OC(1) .CRCI) ,OKG(l> .TCk(l) .TbKtttU ,P<I> 
c  OC  IS  THE  BEAM  CURRENT  IN  UNITS  UF  100  MICHU  AMPS 

C  CR  IS  THE  COUNT  RATE  IN  PHOTONS  PEH  StCONU 

C  BKG  IS  THE  BACKGROUND  COUNT  RATE  IN  PHOTONS  PER  SECOND 

C  TCR  IS  THE  SIGNAL  COUNTING  TIME  IN  SECUNUS 

C  TbKG  IS  The  BACKGROUND  CO  UTINb  TIME  IN  SECONDS  (DARK  COUNT) 

C  P  IS  iHE  PRESSURE  IN  MILLI TOHR 

c  calculation  of  avehage  intensity  and  standard  deviation  at 

c  EACH  ENERGY 

SCRSO.O 
DSCR*O.U 
DU  25 ' lai ,NI 

SLRs&CR* (CR ( 1 ) -OKG (1) )/(BC(I)*P(l) ) 

25  DSCR*USCRi ( CR ( 1 ) / TCR ( I ) *OKG ( I ) / TbKG ( 1 > ) / (BC ( I ) «bC ( I )  *P  1 1  >  *P  (1  > ) 
S(U)*5CR/N1 

C  S(J)  IS  The  AVERAGE  C  ONT  HAIE  AT  THE  ENERGY  C0RRESP0N1NG  TO  J 

D(D)*USURT(USCR)/Nl 

C  0(J)  IS  The  stand  AD  DEVIATION  ASSOCIATED  wITH  S(J) 
c  calculation  of  the  statistical  weighis 

IF (JG.  tU.2)  GO  TO  30 
2b  S ( J) aS  ( J) *0E ( J) 

c  S(vl>  is  the  p  RUUCT  OF  THE  AVERAGE  count  hate  and  the 
c  energy  corresponding  to  u 

D l J) °D ( J) °BE ( J) 

30  CUNTInuE 
KT  a  2 

IF(KT.EU.i)  GO  TO  36 

31  AVSUaO,U 

DU  32  J=1.NJ 

32  AVSD»AvSD*D(U) 

AVSD®AVSD/NJ 

C  A VSD  IS  USED  TO  CALCULATE  THE  WEIGHTS  AND  IS  THE  AVERAGE  OF 

C  The  SlAhUAKU  DEVIATIONS  FUR  THIS  SET  OF  DATA 

DU  35  Jai,NU 
35  wl (J)=(AVSD/D(JI >**2 
GO  TU  4U 

3b  DU  37  Jai ,NJ 
37  WT(J)  a  1.0 
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T able  4.  Continued 


„.0066_ 

0067 

0068 

0069 

0070 

.  00 7J _ 

00  72 
0073 


0074 

0075 

0076 

--0077 

0078 

_M2S _ 

0080 

0061 

0082 

0083 

0084 

__0ft85 _ _ 

0086 
0087 
0088 
0089  . 

0090 

_ 009.1 _ 

0092 

0093 

0094 

0095 


0096 

0097 


0098 


0099 

0100 

.  0101  _ 

0102 

0103 

0104 

0105 

0106 

_0107 


C  WT  ( J)  15  IE  STAl  I5T  1LAL  *E10MT  ASSOCIATED  WITH  S(J) 

*0  CONUNUt 

C  CALCULATION  OF  ABSCISSA  VALUE5 

UU  55  J*  1 1 Nj 

if  ( jo.tu.n  uo  ru  bo 
imjo.eu.2)  oo  ro  51 

50  FOUt(J)3ALUU  (bt(JI) 

bO  TO  55  _  _ 

51  FUUE(O)>1.0/Ut<0) 

55  CONTINUE 

C  FUbE(J)  IS  THE  ABSCISSA  USED  F OH  THESE  SETS  OF  DATA 

C  NtlbHTEO  LtAST  50UAHE5  FIT 

C  CALCULATION  OF  SLOPE  AND  INTERCEPT  ANU  RESPECTIVE  STANOARO 

C  UEV1AT 1UN5. .  _ _ _ _ _  .  ... 

SwT=U.O 

SwTXBV.O 

SWTXAB0.0 

5wTV*0,0 

SWTYYB0.0 

.  S»TAY»0tO  __  . . . . . 

SmwT*u,o 
SWWTXb0.O 
SwwTAAsO.U 
UO  60  J»1*NJ 
SmTbSmT*»T (0) 

.  ..  SWTX«5wTA*wT(j)*fUbEU) _  _ 

SwTXAb5wTXA*WT(J) *F08t ( J) *F06t ( J) 

5bTYbSwTY*wT (0)*S(J) 

SwTYYBSwTYY*wT (J)*STJ)«S(J) 

54TXY»S4TAT*WTUI*StJI*F0BE<Jl 
SwwTB5WWl *WT ( JJ *WT IJ) 

_  SwwTX«SwwTA*wT(J>*WTtJ)*F08e<J> 

60  S«wTAAbS«W1XA»wT ( J) *WT ( J) *FOtit ( J) •F06t (J) 
UEN(R>*SbT»SwIAX-SwTX*SwTX 
SLP(R>b«SbT*Sw1XY-SwTX*SbTY>/UENTK) 

C  SLP(K)  is  THE  SLUPE  OF  THE  OAT A  SET  R 
CtP  (R)=ISwTY#SbT  AX**SwT  X*S*TXY )  /OEN(K) 

_C _  CtPTRI.  IS  THE  INTtHCEPT  OF  THE  OATA  StT  R  _ 

C  CALCULATION  OF  STANDARD  DEVIATION  FOH  SLOPE 

SY (R) “DSUKT { (SwT YY- (SWTY*SWTY*SwTXX-2iU*SWTXY#SwTX,SWT Y 
l*SwT*5WTAY»SWTXY>/UEN<R) )/(N0-2) > 

USLPtRJ  »SY (K) *0SUHT <SWT«SWT*SWWTXX*S*TA»SWTX*SWWT 
1-2«0#5WT#SWTX*SWWTXJ /OEN(K) 

c  __  DSLP(R)  l5  THE  STANDARD  DEVIATION  ASSOCIATED  WITH  SLP(R)_  ... 
DEEP (R) «SY (R) *USUHT l$WTXBSWrA«SWWTXX«SWTXX*5WTXX*SWWW.O* 

1  SWTX»SWTAX*SWWTXJ/0EN(R) 

C  DCEP(R)  IS  THE  STANUARU  DEVIATION  ASSUCIAIEU  WITH  CEP(R) 

If (J6.EU.1)  60  TO  92 

91  HLFE0B1 aO/tO 

OU  TO  93  .  _ _ _ _ 

92  CONTINUE 
REFE0BALOb  <t0) 

93  RNORMtK)  •  SIOMAO/ tSLPTK) *REFE0*CEP IR)  T 

00  94  0B1 t NO 

HETHETOT  ■  HNORH(K)  •  S(J) 

94  UbETHt ( J) BdETHE (O) *USUHT ( ( (OSLP (R) *USLP (R) *  HEFEO  •  REFEO 


•OCE 
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0100 

0109 

0110 


0111 

0112 

0113 

0114 

_U115 _ 

0116 

0117 

one 

0119 

.. 

0121 

0122 

0123 

0124 

012b 

0126 

0127 

0126 

0129 

0130 

0131 

0132 

0133 

0134 

013b 

0136 

0137  . 

0136 

0139 

0140 

0141 

0142 

0143 


Table  4.  Continued 


1P(K)»UCEP<K>  )*HNUHM<K)*HNORM(X)/USLP(XT*  HEFEO  *CEP  <  K )  )  **2)  )  *0  ( 

2J)*DTJ)/<S<J)*S<J>>) _ , _ ' _  ...  .... 

C  BtTHtlj)  lb  Tht  ABSOLUTE  EXCITATION  CROSS  SECTION  FOR  A 
C  RECIPROCAL  ENtHliT  PLOJ 

C  BETHt (J)  IS  THE  ORDINATE  OF  A  bEIHE  PLOT  WHEN  J<>»1  AND 
C  IS  EUUAL  TO  SltiMA*E/U*PI*A*A*R) 

C  RNORM  IS  A  CONSTANT  MHiCH  PUTS  THE  UA  f A  ON  AN  ABSOLUTE  SCALE  BASED 

_C _  ON  THE  ORDINATE. SIOMAO  AT.  tNEK6T_iO. . .  .. 

C  ■  UBETHE  IS  THE  STANDARD  DEVIATION  ASSOCIATED  WITH  THE  value 
C  OF  BETHETJJ 

199  PRINT  1000 

DU  200  JB1 tNJ 

200  PRINT  400(BE<J) fFUBEIJT  *S(J) *U(J) tWT(J) tBETHE(J) iDBETHE(J) 

c  IF  nj  <  IKHAXI  the  HEMAININO  .VALUES  OF  BEf  BETHEi  AND  UBETHE  AHE_ 

C  SEI  EUUAL  TO  ZERO 

IRMAX2  b  IhHAX  •  1 

JO  «  NJ  ♦  1 

DO  16  J  =  J0» 1RMAX2 

BE ( J)  s  0,0 

_  .  BtJHETj)  ■  0*0  _ _  _ 

16  DBETHE(J)  »  0.0 

C  NOW  «E  CHANGE  notation  to  utilize  a  I  mu  dimensional  array 

00  17  J  *  1*1 RMAX 
BEX (JtK)  B  BE ( J) 

BETHX ( Jt  X )  «  UEThE(J) 

17  DBETMXIJiXJ  ■  UBtTME(J) _  .  _  .  .  ... 

6b  CONTINUE 

PRINT  1010 
DU  201  K-ltNK 

201  PRINT  401 tbLP(K)  t DSLP (X) tCEP(K) tDCEP (X) 

PRINT  60b tEO  tSXGMAO 

206  Print  ioio 

DU  202  KB1,NX 

SLP (K ) sRNORM (X) *SLP (K ) 

DSLPTX) =RNORM  <  X ) *USLP ( X ) 

LEP (K  > sRNOHM (K ) *LEP(K) 

ULEP(X)aHNORM(K)*UCtPIX) 

202  PRINT  401tSLPlK)»USLPU),CEP«MT  tUCEPiXJ 
20S  CONTINUE 

C  CHANGE  ENERGY  MATRIX  TO  VECTOR  AND  SAVE  ARRAY  POSITION 
JK  »  1 

UO  301  J  b  i,ikmaa 
DO  301  X  a  UlCMAX 

XBETJX)  =  bEX( J»x)  .  _ 

XBETH ( JK)  a  BETHX (JtK) 

XDBETHIJM  >  DBElHXlJtX) 

IRUM(JK)  ■  J 
1LOLTJK)  »  K 
JX  a  JK  ♦  1 
301  CONTINUE 

C  XBE(JX)  Is  THE  ENERGY  MATRIX  IN  VECTOR  FORM 

C  XBETH(JX)  lb  THE  UKUINATE  MATRIX  IN  VECTOR  FORM 
C  XDBETH(JX)  IS  THE  bTD.  OEV.  MATRIX  In  VECTOR  FORM 

C  IHOM(JK)  IS  THE  HUM  POSITION  VALUE 

C  ILULTJK)  IS  THE  COLUMN  POSITION  VALUE 

C  NO*  ARRANGE  ARRAY  FROM  SMALL  TU  LARGE 
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Table  4.  Continued 


014* 

0145 

_ 01.44 _ 

0147 

0148 

0149 

0150 

0151 

_ftli?i _ 

0153 

—015.4 _ 

0155 
.  0156 
0157 
0158 
0159 

_ . 0160 _ 

0161 

0162 _ 

0163 

0184 _ 

0165 


0166 


0167 

0166 _ 

0169 

0170 

_fll 71 _ 

0172 


0173 
..  0174-.. 

0175 
0176. 
0177 
0176 
0179 
0160 
0161 
— ftl.»2_ 
0183 
ftlM, 
0165 
_0166 
0187 
-111  86 
0169 
Pl»fl 


C  8b AM  ENERGIES  A  HE  ORUEHEU  SUCH  THAT  X6fc  GIVES  THE  6EAM 

C  EHtRGY  AT  COQNUINAItS.  lfiiW..4NU  .l£lik_lliatR£-.ICUU.  CORHESPMUS 

C  TO  THE  ORIGINAL  OATA  SET  NUMBER  AMU  IRUI*  CORRESPONDS  TO  THE 

C  ENtRbY  OAT  A  POINTS  IN  THAI  SET  . 

C  NVAL  IS  THE  NUMBER  OF  VALUES  IN  THE  ARRAY 
NVAL2  «  NVAL  -  1 
DU  302  JK  «  1 iNVAL2 

_ JO.=  OR  *._1  _ _ 

DU  303  A  a  JJtNVAL 

IF  (XBE<UK).Lt,XBE<Kn  .  _.60._T.U.  303 _ _ _ _  _ 

A  ■  XHEWl) 

XdEIK)»X8tlJA»  .  _ _ _ _  _ _ _ _ 

X8E  IJM»A 

_ A  a  XtfETHIM _ 

XbETH(K)  ■  X6ETHUK) 

..  .  XttETIil JK)  »  A  _ _ _ _ _ _ _ 

A  »  XU8ETHIK) 

.  XWBETH(R)  ■  XUBETHI JK) _ 

XUBETHIJR)  a  A 

_ L?iiiy!Li&i _ 

lHOb IK) alROH ( JK) 

_ lHUlHJK)«l  _ _ _ 

I»1C0LIKI 

_ _  „1C0UK)*iq0LTjk)  _ 

ItOLIJK)»l 

303  CONTINUE _ 

302  CONTINUE 

..  _C__— .  NOV  COUNT  NUM6ER._0f_  UKE.EiltRUIti _ _ _ _ 

N  a  0 

..  C  N  COUNTS  LIKE  VALUES -OF.  B£A«_feNEIi6Y»_xee.  .  _ 

DO  304  JK  a  i,NVAL 

_ j£  tJK.EU.l)  GO  TO  320 _ 

IF  (X6E(JK>.EU.CAB£(N)>  GO  TO  304 

_ _C  .  .  IF  XBEU)  £UUALSCABElN)_i_iTlHASALHEAUY  BEEN  COUNTED..  . 

320  N  a  N  ♦  1 

CXbE(N)  ■  XGE(JK)  .  ...  , _  _  .  . 

1LOUNT IN)  ■  I 

C  ICOUNt  IS  IHE  NUMBER  OF  TIMES  EACH  ENERGY t  CXBEt  APPEAHS _ 

K  a  JK  ♦  1 

.  ..  DO.  305  L  ?  KjNVAL  ...  .  ..  . .  . .  . 

IF  ICXBEINl.NE.XbbtL)')  GO  TO  349 

ICOUNI(N)  «  ICOUNT  (N»  .♦_! _ . . 

305  CONTINUE 

349  CONTINUE _ 

304  CONTINUE 

_ .ZENO  »  .  .  _ _  _  _ _ 

LX  *  i 

_  .  IF  (CX6E ILX1  tNEtlEHO) GO.JO.29i . 

1  b  ICOUNT (1)  ♦  1 

: _ DO..  292__OfiJL_»_lt«VAL _ 

JK  a  JKK-  I  ♦  1 

.  .  XbETHUK)  »  XbETH(UKK)  .  ....  ...  .  . 

292  XU6E1H(JM  a  xUbbTHl JKK) 

NVAL  »  NVAL  -  I  ♦ .  1  ..  . .  . 

DO  290  KK  a  2»N 

_  L  ■  KK  -  1 _ : _ 
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Table  4.  Continued 


0191  ICOUNT (L>  ■  1  COUNT  IKK) 

U192 _ 2301  ..CABElL)  g_CABLiem _ 


0193 

N  a  N  -  1 

0194 

291 

CONTINUE  _  .  ...  .  .  _ 

C 

NUM  AVERAGE  BETHE  AND  COMBINE  STANDARO  DEVIATIONS  AT 

C 

EACH  ENERGY  _ _ 

0195 

JJ  ■  I 

0197  LL  «  ICOUNT  ID  ♦  JJ  -  1 

— oivo  sun  «o.o 


0199  SSUM  ■  0.0 

.  . 020.0  00  30/  OK  ■  JJ'LL  .  .  .. 

0201  SUM  ■  SUM  ♦  XBETH  I JK) 

_ MI2 _ S>UM_ ■  SSUM  ♦  XDBETH ( JK1  axDBETH ( JK) 

0203  307  CONTINUE 

0204  ABETHE(L)  «  SUM/1COUNML)  _ 

0205  AUBETH(L)  aOSURT (SSUM) /ICOUNT (L) 

.  0200  .  JJ  ■  LL  ♦  l  _  _ 

0207  306  CONTINUE 


—0208 
0209 
-0210 
0211 
_0212_. 
0213 
0214 
021S 
-.0210 
0217 
0  218 
0219 
0220 
0221 


0223 
.  0224 


0225 

022fc_ 

0227 

0228 

0229 

0230. 

0231 

0232 

0233 


0234 
0235 
...  023fe. 


X _ AUKTHt.-lS  THE  AVERAGE  ORDINATE  KITH  STAND AMD  DEVIATION  ADBETHE 

C  FOR  EACH  DISCRETE  BEAM  ENERGY  CXBE 

C  WEIGHTED  LEAST  SUUAHES  FIT 

C  CALCULATION  OF  ABSCISSA  VALUES 

1)0  309  L*ltN  . 

BE(L)»CXBE(L) 

_ SAL)  ■A.BEIHE  (L) _ 

309  U(L) s*DBETH  (L) 

IF (J6.EU.2)  GO  TO  312 
PRINT  603 

312  CONTINUE 
PRINT  604 

_ DO  316  L!1*N  _ _ _ _ 

IF (JG.EU.l)  GO  TO  315 

310  FOBElL)»i.O/BE(U 
GO  TO  316 

315  FUBE(L)>ALOG  (BEtD) 

316  CONTINUE 

IF  (UG.NEf 3)  GO  TO  .313 _ 

PNINT  603 

313  CONTINUE 

C  AVSD  IS  USEO  IN  THE  WEIGHT  CALCULATION  ANO  IS  THE  AVERAGE 

C  OF  THE  STANUAKU  OEV1AT IONS  FUN  THIS  SE1  OF  DATA 

AVSDsO. 0 

.  DO  306  L»ltN  _ _ _ 

308  AVSD  «  AVSO  ♦  O(L) 

AVSD°AVSD/N 
DO  311  L®1 iN 

311  WT  (L>  a  (AVSI)/0(L)  I **2 
PRINT  1020 

_  00  351  L  a  1 ,N  _ 

351  PH1NT  S04tbE(L) tFOBE (L) *S(L) »0(L) *WT(L) 

C  CALCULATION  OF  SLOPE  AND  INTERCEPT  ANO  RESPECTIVE  STANDARD 
C  DEVIATIONS 

SwT»0.0 
SaTX«U.U 
SwTXA»0.o 
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Table  4.  Concluded 


0237  SKTY*0.,0 

0238 _ S"TYYB0.0 

0239  SMTXr*0.0 


024a. 

0241 

0242 

0243 

0244 _ 

0245 

0246  .  . 

0247 

0246. 

0249 

_ 0250 _ 

0251 

_  0252. _ 

0253 

0254 _ 

0255 

.  0256 _ 

..  0257. _ 


0256 


SMHTbU.O 

S«IWTXb0.« 

5wbTAXb0.0 
DU  350  LBltN 

_ .Sd.?.&.8l*6T  ID _ 

S«TXBSvT X*nT CL» *7 Obfc <L) 

ShTXXbSMTXX*WT  (Li  *FUBE  (L)'*F0Bt  (L)  . 
SmTYb5WTY«wT(U*S(U 

SNTYYBSBlYY«BT(U*S(L)*S<U  -  _  -  - 

SwTXYBS8l XY+WT(L)*S(L) *F0BE (L) 

_ &MaI?5MkT*llT  1LL*hT.(U  _____ - 

S6wtxbsmmtx«mt  <u  *«r  tu*Fout  iu 
350  SHWTXXbSNWTAX*HT  (L)  *VT  (L)  *FOBE(LJ  4FQBEIL)  ...  . 
DEmbSnT*ShTXX-Sm  T  X*S*T  X 

_  .  SlMB<SWT«Sl»TXY-SwTX*&Bm/DEM  _ 

C  SLM  15  THE  SLOPE  OF  THE  DATA 

_ CtK^5ttLY*J.«TXX-58lAlSMimZUtM - 

C  CEH  15  THE  INTERCEPT  IF  THE  DATA 

_  .  SM«U5WHT ( (5hTYY“ tSwTY*5WTY*SW rxX”2.0*5»TXY*SWTX*SBTY 

1  •SRT*SbTXY*SWTXY>/DEN)/(N-2)> 
DSLMB5M*05URTISWT*5WT*5B*TXX*S*TX*SBTX*SBBT 
1  -2.U*S«*T*SBTX*SWbTA»/DEH 

JC _ ttSLM-IS  THE  STANDARD  DEVIATION  ASSOCIATED  WITH  SLM 

OC£MbSM*USuRT  <SwTX*SwTA*SRBTXX*S«TXA*SVTXX*SWbT-2.04 
.  .  1  5MTA*5HTXX*SHUlX)/ULli 

C  DCEM  IS  THE  STANDARD  DEVIATION  ASSOCIATED  «1TH  CEH 


0259  _ _ _  PRINT  550 •  SLMiUSLM 

0260  PRINT  551  *  CEH.ULEM 

0261 _ IfJue.EUflJ  .  GUJLSLJ53 _ 

0262  IF(JG.GT.2)  60  TO  353 

_ 5263 _  CONST  «  4.788  t-lb 

C  CONST  IS  EUUAL  TO  4*PI*A*A*R  IN  UNITS  OF  EV*CM#CM,  WHERE 

_  _ t  A  IS  THE  FIRST  BOHR  RADIUS  AND  R  15  THE  RYDBERG  ENERGY 

0264  DO  354  L  ■  1  *N 

0265 _ FOBE(L)  B  ALUtolBElEU _ 

0266  S(L)  *  S(U*BE(L) /CONST 

0267 _ ..354  D<U  ■  U  (U  •6E.IL)  /CONST  .  _  ... 

0268  JG  ■  JO  ♦  1 

0269  .  GO  TO  ' 316  .  . . . 

0270  353  STOP 

0271 _ END _ 
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By  following  the  derivation  of  Beers  (Ref.  15)  for  a  least-squares 
fit  and  noting  the  derivation  of  Bacon  (Ref.  16)  for  the  weighted  slope 
and  intercept,  equations  for  the  weighted  standard  deviations  of  the 
weighted  slope  and  intercept  were  derived,  the  details  of  which  are  listed 
in  the  appendix  of  Ref.  17.  Beginning  with  the  standard  deviation  of  a 
random  count  as  the  square  root  of  the  count,  the  computer  calculation 
computes  the  statistical  weight  associated  with  each  energy  and  performs 
the  weighted  least-squares  fit  for  each  energy  scan.  The  coordinates 
are  a  versus  E-1  or  aElAira^R  versus  in  E  as  desired.  Next,  a  normali¬ 
zation  constant  is  computed  such  that  the  straight  line  determined  by  the 
weighted  least-squares  fit  passes  through  the  ordinate  value  used  as  in¬ 
put  data  at  10  keV,  and  the  observed  intensity  measurements  determined 
as  a  function  of  electron  beam  energy  are  put  on  an  absolute  scale  by 
multiplying  each  value  by  this  normalization  constant.  After  each  set 
of  data  is  thus  scaled  to  give  absolute  apparent  cross  sections,  or 
(jEM^a^R  for  a  Bethe  plot,  the  data  are  sorted  to  give  all  cross-section 
values,  or  ctE / for  a  Bethe  plot,  and  their  respective  standard 
deviations  occurring  at  each  electron  beam  energy.  After  averaging  the 
ordinate  values,  a  weighted  least-squares  fit  is  performed,  and  the  ordi¬ 
nate  and  abscissa  values  as  well  as  the  slope  and  intercept  are  printed 
along  with  their  respective  standard  deviations  and  statistical  weights. 

Thus,  the  observed  intensity  measurements  for  each  spectral  line  deter¬ 
mined  as  a  function  of  electron  beam  energy  were  used  with  the  10-keV 
data  of  the  Nil  lines  and  the  Njj^l-MO,  0)  band  to  obtain  absolute  values  of 
the  apparent  cross  sections.  The  values  thus  obtained  for  the  apparent 
cross  sections  of  the  Nil  3F°  and  Nil  states  produced  by  dissociative 
excitation  of  N2  and  NO  are  listed  in  Tables  5  and  6,  respectively. 

1 

The  partial  Grotrian  diagram  shown  in  Fig.  22  indicates  qualitatively 
the  possible  cascade  contributions  to  the  Nil  3F°  and  Nil  1D  levels  (Ref.  18). 
The  possible  existence  of  these  transitions  induced  by  electron  beam  exci¬ 
tation  is  at  present  unknown,  and  thus  no  cascade  corrections,  if  they  do 
indeed  exist,  can  be  made  on  the  present  cross-section  values. 

The  apparent  cross  sections  of  the  Nil  level  produced  by  dissoci¬ 
ative  excitation  of  N2  and  NO  result  directly  from  the  optical  cross-section 
measurements  since  the  Nil  level  radiates  by  only  the  observed  tran¬ 
sition  (Fig.  22).  The  Nil  ^F°  level,  however,  has  three  competing  radi¬ 
ative  transitions  (Fig.  22),  and  the  optical  cross  section  at  10  keV  was 
divided  by  the  branching  factor  of  0.  155  as  determined  from  Ref.  18  in 
order  to  obtain  the  apparent  cross  section.  Consequently,  the  apparent 
cross  sections  possess  the  accuracy  of  the  life  time  calculations  of  Ref.  18. 
Since  the  life  time  is  the  same  for  the  NI  3F°  level,  regardless  of  molecular 
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Table  5.  Energy  Dependence  of  the  NII(3F°)  Apparent  Cross  Section 


n2 

NO 

E(keV) 

g(10"^  cm2) 

oElO3 
4Tia02  R 

E(keV) 

o(10~^  cm2) 

OE103 

4Tra02R 

3.6 

isa 

4.0 

2.8 

2.3 

2.8 

HI 

5.0 

2.3 

2.4 

2.4 

3.0 

6.0 

1.9 

2.4 

7.0 

2.0 

3.0 

7.0 

1.7 

2.5 

8.0 

1.8 

3.0 

8.0 

1.5 

2.6 

10.0 

1.5 

3.0 

10.0 

1.3 

2.7 

12.0 

1.2 

3.0 

12.0 

1.1 

2.7 

15.0 

0.92 

2.9 

15.0 

0.89 

2.8 

18.0 

0.89 

3.4 

18.0 

0.75 

2.8 

20.0 

0.89 

3.7 

20.0 

0.68 

2.8 

Table  6.  Energy  Dependence  of  the  NIM1  D)  Apparent  Cross  Section 
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specie  involved,  the  ratjo  of  the  NIKN2;  5016.4  A)  optical  cross  section 
and  the  NII(NO;  5016.4  A)  optical  cross  section  is  identical  to  that  of  the 
apparent  cross-section  ratio,  and  the  dependence  on  the  accuracy  of  the 
life  time  calculation  is  eliminated.  The  ratios  thus  obtained  are 


ff[NII(N2;5016.4  A)]/ff[NII(.NO;  5016.4  A)]  =  1.16 

and 

a[NII(N2;  3995.0  A)] /o{NII(N0;  3995.0  A)]  =  1.22 
determined  at  an  electron  beam  energy  of  10  keV. 

I 

By  using  the  Bethe-Born  analysis  of  Section  1.2,  the  values  of 
several  parameters  can  be  derived.  For  optically  allowed  transitions, 
the  slope  of  a  Bethe  plot  is  equal  to  the  dipole -matrix- element  squared, 
for  the  transition;  the  constant  in  Cn  is  equal  to  the  ratio  of  the 
intercept  to  the  slope;  the  constant,  Ka0,  is  derived  from  the  value  of 
Cn  using  Eq.  (30);  and  the  optical  oscillator  strength,  is  defined  by 
Eq.  (12).  First-order  forbidden  transitions  yield  the  constant  bn  by 
dividing  the  slope  of  a  reciprocal  energy  plot  by  47ra2R  which  is  equiva¬ 
lent  to  using  the  expression  bn  =  aE/ 47raQR.  For  the  N2(l-)(0,0),  the 


53 


AEDC-TR -75-77 


following  parameters  were  calculated: 

VI2  =  0.353  ±  0.009 

n 

C  =  0.035  ±  0.009 

n 

Kao  =  0.110  ±  0.001 


and 


f  =  0.49  +  0.01 

n 


4.0  DISCUSSION  AND  CONCLUSIONS 


The  dependence  of  the  Nr^l-KO.O)  band  emission  cross  section  on 
the  energy  of  the  incident  electrons  was  in  excellent  agreement  with  the 
results  of  Ref.  13,  but  only  after  steps  were  taken  to  eliminate  secondary 
electron  excitations.  At  low  electron  beam  energies,  secondary  excita¬ 
tion  processes  due  to  backscattered  electrons  from  the  Faraday  cup  were 
appreciable  and  contributed  as  much  as  13  percent  to  the  observed  inten¬ 
sity  at  4  keV  but  were  negligible  near  20  keV.  The  30 -v  negatively 
biased  grid  on  the  Faraday  was  more  than  sufficient  to  suppress  these 
secondary  electrons. 

The  theoretical  calculation  of  the  N2(l-)(0,0)  band  profile  described 
in  Section  3.  0  was  in  excellent  agreement  with  the  experimentally  ob¬ 
tained  profile.  The  ratio  of  the  R  branch  peak  height  to  the  P  branch 
peak  height  was  0.  305  for  both  the  theoretically  and  experimentally  de¬ 
termined  profiles.  From  this  measurement  and  those  of  Ref.  5  at  other 
temperatures,  it  is  believed  that  the  procedure  used  in  this  paper  is 
accurate  to  within  less  than  0.  5  percent. 


The  apparent  cross  sections  of  the  Nil  °F°  state  produced  by  the  dis¬ 
sociative  excitation  of  N2  (Fig.  23)  and  NO  (Fig.  24)  over  the  electron 
beam  energy  range  of  4  to  20  keV  and  plotted  against  the  reciprocal 
of  the  electron  beam  energy  show,  as  do  the  Bethe  plots  (Figs.  25  and 
26),  that  the  transition  is  consistent  with  the  Bethe-Born  equation  for 
forbidden  transitions.  The  apparent  cross  section  of  the  Nil  ■'•D  state 
produced  by  the  dissociative  excitation  of  N2  (Figs.  27  and  28)  is  also 
consistent  with  the  Bethe-Born  approximation  for  a  forbidden  transition. 
The  apparent  cross  section  of  the  Nil  1D  differs  significantly  from  the 
Bethe-Born  equation  below  about  8  keV  as  shown  in  Figs.  29  and  30. 
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Figure  23.  Energy  dependence  of  the  molecular  nitrogen  NII(3F°) 
apparent  cross  section. 
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Figure  24.  Energy  dependence  of  the  nitric  oxide  NII(3F°) 
apparent  cross  section. 
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Figure  25.  Bethe  plot  of  the  molecular  nitrogen  Nll{3  F°) 
apparent  cross  section. 
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Figure  26.  Bethe  plot  of  the  nitric  oxide  NII(3F°) 
apparent  cross  section. 
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Figure  27.  Energy  dependence  of  the  molecular  nitrogen  NIICD) 
apparent  cross  section. 
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Figure  28.  Bethe  plot  of  the  molecular  nitrogen  NIK1  D) 


apparent  cross  section. 
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Figure  29.  Energy  dependence  of  the  nitric  oxide  NIK1  D) 
apparent  cross  section. 
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Figure  30.  Bethe  plot  of  the  nitric  oxide  Nl  I C  D) 
apparent  cross  section. 
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The  slopes  as  determined  by  the  Bethe  plots  are  [6. 1  ±  8.  3]  x  10“*\ 

[3.5  ±  0.  5]  x  10-4,  o[2. 4  ±2.5]  x  10~4&  and  [-6.7  ±  7.0]  x  10‘5  for 
the  NII(N2;  5016.4  A),  NII(NO;  5016.4  A),  NII(N2;  3995.0  A),  and 
NII(NO;  3995.  0  A)  transitions,  respectively.  These  values  along  with 
the  respective  standard  deviations  show  a  near  zero  or  zero  slope  indi¬ 
cating  that  the  Bethe -Born  equation  for  a  forbidden  transition  is  appli¬ 
cable  . 

The  nonlinearity  below  8  keV  on  the  NII(NO;  3995.  0  A)  Bethe  plot 
(Fig.  30)  may  be  due  to  the  existence  of  two  competing  processes  pro¬ 
ducing  the  same  radiating  Nil  ■’-D  state.  First,  consider  the  possibility 
that  the  Nil  state  is  produced  by  two  competing  dissociative  exci¬ 
tation  processes  where  one  process  is  optically  allowed  and  the  other 
forbidden,  the  optically  allowed  process  being  dominant  at  the  lower 
energies.  Then  as  the  electron  beam  energy  is  increased,  the  cross 
section  of  the  optically  allowed  process  becomes  sufficiently  small  to 
enable  the  forbidden  process  to  dominate,  thus  producing  the  observed 
energy  dependence.  Next,  consider  the  situation  where  another  excited 
state  (Fig.  22)  as  well  as  the  Nil  is  produced  by  the  dissociative 
excitation  process,  the  other  excited  state  subsequently  cascading  into 
the  Nil  level.  If  the  higher  excited  state  is  produced  by  an  allowed 
process  and  the  Nil  state  by  a  forbidden  process,  the  observed  be¬ 
havior  could  result.  Thirdly,  it  is  possible  that  the  Nil  1D  state  is  not 
excited  by  the  incident  electrons  but  is  due  only  to  two  cascading  tran¬ 
sitions,  one  of  which  is  produced  by  an  allowed  and  the  other  a  forbidden 
dissociative  excitation  process.  Again  the  observe  energy  dependence 
could  result.  Note,  however,  that,  in  all  three  cases,  there  is  one 
allowed  and  one  forbidden  process.  For  the  cases  where  both  processes 
are  forbidden,  only  a  different  value  of  bn  would  result.  Little  can  be 
said  about  the  validity  of  the  cascading  arguments,  however,  as  the 
population  of  these  states,  which  may  cascade  with  the  Nil  state,  by 
electron  collision  is  at  present  unknown. 

Potential  energy  curves  for  the  electronic  states  of  N2  and  NO  exist 
(Ref.  19)  only  for  the  lowest  electronic  levels  and  are  not  adequate  to 
provide  information  regarding  the  excitation  channels  which  produce 
Nil  and  Nil  ^F°  excitations.  In  principle,  however,  it  is  well  known 
that  when  a  molecule  in  its  ground  electronic  state  undergoes  a  collision 
with  a  high  energy  electron,  excited  electronic  states  may  be  produced, 
some  of  which  may  lead  to  dissociation  of  the  molecule,  and  the  excitation 
processes  obey  the  Franck-Condon  principle  as  described  in  Section  1.0. 
Therefore,  excitations  may  occur  to  stable  electronic  states,  to  elec¬ 
tronic  states  from  which  the  molecule  will  dissociate,  and  to  electronic 
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states  where  both  stable  molecules  and  dissociation  results.  A  number 
of  electronic  states  of  the  last  two  types  may  exist  and  may  result  in 
the  Nil  excitations  of  interest.  Thus,  accurate  potential  energy  curves 
are  needed  to  aid  in  the  understanding  of  the  physical  processes  in¬ 
volved  in  these  excitations. 


Finally,  as  an  application  of  the  results  of  this  study  we  shall  con¬ 
sider  the  use  of  the  electron  beam  fluorescence  technique  for  specie 
density  measurements  of  N2/He  and  NO/ He  binary  mixtures.  Con¬ 
sidering  the  excitation  of  He  by  electron  collision,  it  is  well  known  that 
the  process 

e"  +  IIe(l  ’S)  *  HeCl’P)  *  e‘ 

L  Hef2*S)  -r  Mom 6  A)  (57) 


is  desirable  for  utilization  for  diagnostics  purposes  because  of  the  ease 
with  which  it  is  excited.  It  is  recalled  that  both  N2  and  NO  produce 
radiation  at  5016  A  from  Nil  radiation  transitions  resulting  from  elec¬ 
tron  dissociative  excitationoprocesses.  Consequently,  the  measured 
radiative  intensity  at  5016  A  observed  from  a  binary  mixture  of  either 
N2/He  or  NO/ He  results  from  both  He  and  either  N2  or  NO  excitations 
and  is  not  a  unique  measure  of  an  individual  specie  density.  The  contri¬ 
bution  to  the  measured  intensity  by  the  individual  species  is,  of  course, 
dependent  not  only  on  the  specie  density  but  also  on  the  individual  exci¬ 
tation  cross  section.  The  following  section  will  consider  such  mixtures 
and  the  magnitudes  of  the  individual  specie  excitation  cross  sections  for 
the  purpose  of  determining  binary  mixture  compositions  for  which  unam¬ 
biguous  and  accurate  density  measurements  are  possible. 


Consider  the  interaction  of  an  electron  beam  of  current  I  and  cross- 
sectional  area  (a)  with  a  gaseous  binary  mixture  consisting  of  specie  1 
and  specie  2  in  their  ground  states  g  such  that  excited  states  i  and  j, 
respectively,  are  produced  and  a  steady -state  photon  emission  rate  is 
observed  at  a  wavelength  The  volume  rates  of  production  of  atoms 
are  (I/a)ni  <rL  for  the  state  i  and  (I/a)n?  ctL  for  the  state  j,  and  the 

S  s1  J  6J  1  1  1  1 

•  .  n  4-  A  n  A-f  M  a]  t  ^  1  t  r  A  J  /%  il  +  4  *-1  -V*  A  rCl  IT  AM  Vvt  T  )  A  M  -  n  A  .  M  . 


volume  rates  of  radiative  depopulation  are  given  by  ?  Ak  ni 


=  Ai 


m  for 


the  state  i  and  £  Af„ 
r  Jr 


nj 


2  2 


1 


=  Ai  n.j  for  the  state  j,  where  Og^  is  the 


J  3 


excitation  cross  section  of  specie  1  for  the  g  -»  i  transition,  n_.  is  the 
specie  1  ground  state  number  density,  and  A^  is  the  specie  1  Einstein 
transition  probability  for  a  spontaneous  i  -  JL  transition.  When  it  is 
assumed  that  no  other  processes  are  involved,  in  particular,  radiative 
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cascading  and  collisional  quenching  effects,  the  steady-state  number 
densities  in  states  i  and  j  are  given  by 


n}  =  (I/a)nJaJ./A  ? 


and 


nf  - 


(58) 


(59) 


For  i  -*  s  and  j  -*  t  transitions,  multiply  each  equation  by  the  appropriate 
Einstein  transition  probability  to  obtain  the  photon  emission  rates  per 
unit  volume. 


51  .  A‘„„!  .  U'ajn'^/s!, 

52  -  Af,»f  - 


(60) 

(61) 


where 


ft?  =  A?  'A.1 

^IS  IS  I 

is  defined  as  the  branching  factor. 

The  observed  photon  emission  rate  is  given  by 

<S‘  *  s2>V.bs«,..„o.  -  (S'  +  SZ)L. 


where  L  is  the  observed  electron  beam  length.  Therefore,  the  observed 
photon  emission  rate  S  is  given  by 

s  -  s'  +  s2  _  ♦  ngffgj  pftj  ■  (62) 

Note  that  neither  the  photon  emission  rate  for  specie  1  nor  specie  2  can 
be  determined  independently.  Therefore,  in  a  binary  mixture  of  N2/He 
or  NO/ He,  for  which  the  He  mole  fraction  is  small,  the  observed  inten¬ 
sity  of  the  He(5016  A^P  "*  2*S  transition,  which  can  easily  be  obscured 
by  the  NII(5016.4  A)*D  ■*  1p°  transition,  cannot  be  distinguished  as  coming 
from  the  He  or  the  Nil  transition.  For  the  cases  of  N2/NO,  N2/N,  and 
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1  2 

NO/N  mixtures,  Eq.  (62)  is  changed  only  in  that  /3^g  =  j3jt,  and  the  sub¬ 
scripts  i  and  j  refer  to  the  same  state  as  do  s  and  t.  Thus,  the  observed 
intensity  from  an  excited  Nil  state  is  indistinguishable  as  coming  from 
N2,  NO,  or  N. 

The  technique  demonstrated  herein  applies  to  flow  diagnostics  as 
well  as  determining  interesting  cross  section  behavior.  The  range  of 
mixtures  of  N2/He  and  NO/He  for  which  one  can  use  the  5016-A  He  line 
for  flow  diagnostics  is  determined  by  the  relative  densities,  transition 
probabilities,  and  cross  sections.  By  using  the  transition  probabilities 
of  Ref.  18,  the  optical  cross  section  for  He  5016  A(3*P  -*  2^S)  extrapo¬ 
lated  to  10  keV  from  Ref.  20,  and  the  present  optical  cross  sections 
for  the  Nil  5016.  4-A  transitions,  the  ratios  of  He  5016-A  intensity  to  the 
total  observed  intensity  at  10  keV  as  a  function  of  n^e / nj^  and  njje/nNC) 
was  calculated  and  is  shown  in  Fig.  31.  It  is  seen  that  unless  the  number 
density  of  helium  is  several  times  that  of  N2  or  NO,  a  substantial  part 
of  the  observed  intensity  is  due  to  the  Nil  transition.  Thus,  previous  and 
future  local  specie  density  measurements  using  the  He  5016-A  line  must 
include  corrections  for  Nil  5016.  4-A  contributions  to  the  observed  inten¬ 
sity. 


nH  e^X 


Figure  31.  Ratio  of  helium  5016  A  to  total  observed  intensity  at 
10  keV  as  a  function  of  number  density  for  N2/He  and 
NO/He  mixtures. 
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NOMENCLATURE 

Einstein  transition  probability  for  spontaneous  i  -*■  s 
transition 

Electron  beam  cross-sectional  area 
Bohr  radius 

Upper  electronic  state  of  nitrogen 

Constant  defined  by  Eq.  (33)  and  appearing  in  the  excitation 
cross-section  relation  for  optically  forbidden  transitions 

Constant  defined  by  Eq.  (30)  and  appearing  in  the  excitation 
cross-section  relation  for  optically  allowed  transitions 

Constants;  1.  1909  x  10~12  watts-cm2/sr  and  1.4380  cm-K, 
respectively 

Speed  of  light 

Differential  cross  section  for  excitation  from  ground  state 
to  state  n  with  incident  particle  scattered  into  solid  angle 
element  du 

Kinetic  energy  of  incident  electrons 
Excitation  energy 

Excitation  energy  from  ground  state  to  state  n 
Electronic  charge  and  electron 
Optical  oscillator  strength 
Generalized  oscillator  strength 
Photon  of  energy  hi/ 

Planck's  constant  divided  by  2 tt 
Initial  momentum  of  the  incident  electron 
Fined  momentum  of  the  incident  electron 
Momentum  transfer  of  the  incident  electron 
Electron  beam  current 
First  ionization  threshhold 

Measured  relative  intensity  of  the  tungsten  ribbon  standard 
lamp  using  the  optical  system  to  be  calibrated 
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Ka0  Constant  defined  by  Eq.  (28) 

L  Observed  electron  beam  length 

M  Reduced  mass  of  the  incident  electron -molecule  system 

_  Z  Z  Z 

M  Vector  with  components  £  ejxj  ,  £  ejyj  ,  and  ^  ejzj 

j=l  j  =  l  ‘  j  =  l 

M„  Dipole -matrix-element  squared  for  transition  from  ground 

state  to  state  n 

m  Electronic  mass 

NI  Neutral  atomic  nitrogen 

Nil  Singly  ionized  atomic  nitrogen 

NjU-)  First  negative  system  of  ionized  molecular  nitrogen 


N2<2+)  Second  positive  system  of  molecular  nitrogen 

NjjU-Hv,  v")  v  -*  v'  vibrational  band  of  NjjU") 

N^  Spectral  radiance,  watts/ cm1 & -sr- cm 


n 

A 

n 

Pnm 

P 

R 

-*■ 

r 


S 

S(A ) 
T 

tb 

V 


xl4 


Number  density 
Unit  vector 

Probability  of  ann->m  transition 
Electron  momentum 
Rydberg  energy 

Position  of  the  incident  electron  relative  to  the  center  of 
the  atom 

Position  of  the  atomic  electrons 
Photon  emission  rate 
Optical  system  relative  sensitivity 
True  temperature 
Brightness  temperature 

Coulomb  interaction  potential  energy  between  the  incident 
and  atomic  electrons 

Velocity  of  incident  electron 

Ground  electronic  state  of  nitrogen 
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yjm 

z 

^is 

e(*,T) 

en 

r}(E') 

®» $ 

BA 

\ 

5n<4?) 

P 

CTij 


^nv 
^nvJM 
V'n^l  •  •  • 

w 


Lower  electronic  state  of  nitrogen 
Component  of  the  position  vector  rj 
Spherical  harmonics 

Number  of  atomic  electrons  in  target  atom 
Branching  factor  for  the  i  -♦  s  transition 
Emissivity 

Atomic  matrix  element  or  form  factor 
Efficiency  or  probability  factor 
Spherical  polar  coordinates 
Spherical  polar  coordinates 

o 

Wavelength,  A 
Electronic  wave  function 
Internudear  distance 

Total  cross  section  for  the  i  -►  j  transition 

Total  cross  section  for  excitation  from  the  ground  state 
to  state  n,  regardless  of  angle  of  scattering  of  incident 
particle 

Transmission  factor  of  standard  lamp  quartz  window 
Vibrational  wave  function 
Total  molecular  wave  function 

,  rz)  Eigenfunctions  of  state  n  and  in  the  coordinates  ?j  of  the 
atomic  electrons 

Solid  angle  subtended  by  scattered  electrons 
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